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region  when  a  high  electric  field  is  applied  parallel  to  the  interface. 

It  is  analogous  in  many  respects  to  the  Gunn  effect,  except  that  this 
mechanism  allows  greater  control  of  device  characteristics.  These 
characteristics  can  be  adjusted  by  varying  the  doping  densities,  the 
layer  thicknesses,  and  the  A1  mole  fraction  in  the  AlGaAs. 

The  mechanism  is  analyzed  using  the  electron  temperature  model 
and  the  Monte  Carlo  simulation.  The  electron  temperature  model  is  exact 
in  the  high  carrier  density  limit,  whereas  the  Monte  Carlo  method  is  valid 
in  the  low  density  limit.  Both  methods  clearly  illustrate  the  degree  of 
control  possible  with  this  mechanism  over  device  characteristics.  Com¬ 
parisons  are  made  between  the  two  models.  Miscellaneous  effects  which  are 
neglected  in  the  models  are  discussed.  These  include  two-dimensional 
effects,  band  bending,  statistical  fluctuation,  and  quantum  mechanical 
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transmission  at  the  interface.  Switching  characteristics  have  been 
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and  memory  devices. 


The  second  portion  of  the  present  work  deals  with  the  band 
structure  dependence  of  high  field  transport  and  impact  ionization  in  GaAs. 
An  understanding  of  impact  ionization  is  important  because  of  its  influence 
on  the  performance  of  avalanche  photodetectors  and  IMPATT  diodes.  It  also 
determines  the  ultimate  performance  limits  of  small  semiconductor  devices, 
such  as  CCDs  and  FETs. 

A  realistic  band  structure  has  been  included  in  a  Monte  Carlo 
simulation  of  high  field  transport  in  GaAs.  The  band  structure  has  been 
calculated  using  the  empirical  pseudopotential  method.  Partly  due  to  the 
lack  of  information  and  partly  for  simplicity,  we  have  made  simplifying 
assumptions  on  the  phonon  scattering  rate,  the  ionization  threshold,  and 
the  ionization  probability.  Unlike  previous  theories  of  impact  ionization, 
the  method  requires,  in  principle,  no  adjustable  parameters  as  long  as  the 
band  structure  and  the  scattering  mechanisms  are  known.  The  calculated 
drift  velocity,  mean  free  path,  and  impact  ionization  rate  are  in  fair 
agreement  with  the  experimental  data.  It  is  found  that  the  contribution 
of  ballistic  electrons  to  the  impact  ionization  rate  is  negligibly  small. 
Within  statistical  uncertainty  we  do  not  observe  the  anisotropy  of  the 
electron  ionization  rate  in  contradiction  to  the  recent  experimental 
results.  Based  on  the  results  of  the  simulation,  a  general  discussion  of 
impact  ionization  is  given. 
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1  .  INTRODUCTION 

The  field  of  hot  electron  physics  is  now  more  than  thirty  years  old. 
Following  earlier  work  on  dielectric  breakdown  in  insulators,  Shockley[l], 
in  1951,  initiated  experimental  investigation  of  hot  carrier  effects  in 
semiconductors.  In  the  60^3  the  field  saw  one  of  its  peaks  with  the 
prediction  and  the  discovery  of  the  Gunn  effect.  The  historical  background 
which  led  to  thi3  discovery  is  discribed  by  Hilsum  in  his  excellent 
paper[2].  The  recent  revival  of  interest  in  this  field  is  due  to  hot 
carrier  effects  which  occur  in  submicron  silicon  MOSFETs  and  GaA3  FETs . 
Oxide  breakdown ,  hot  electron  emission  from  silicon  into  the  gate- oxide, 
velocity  overshoot,  and  ballistic  transport  in  GaAs  FETs  are  examples  of 
such  effects. 

Among  all  the  semiconductors  the  III-'/  compounds  rank  only  behind 
silicon  in  technological  importance,  and  are  second  to  none  with  respect  to 
basic  research  interests.  One  of  the  unique  properties  of  the  III-'/ 
semiconductors  is  their  band  structure  which  makes  possible  the 
construction  of  transferred  electron  devices  (TEDs)[3j-  The  fascinating 
idea  of  electrons  transfering  in  momentum  space  from  one  valley  to  another 
has  attracted  many  physicists  and  engineers[2 ] .  Another  interesting 
property  of  III-V  semiconductors  is  that  they  can  be  epitaxially  grown 
lattice-matched  to  each  other  when  proper  binary,  ternary  o”  quaternary 
alloys  are  employed.  This  property  has  been  fully  utilized  in 
optoelectronics  devices.  From  the  viewpoint  of  carrier  transport, 
heterostructures  can  orovide  an  entire  new  class  of  devices  with  unusual 
transport  properties.  The  concept  of  a  "superlattice"  first  proposed  by 
Esaki  and  Tsui"!]  in  1969  is  one  such  example.  It  is  expected  that  in  the 


•.Uuure  more  and  more  heterostructures  ’//ill  be  utilized  in  various 
semiconductor  devices. 

It  is  possible,  for  example,  to  realize  electron  transfer  in  real 
space  utilizing  heterostructures  as  proposed  by  Hess  and  co-workers] 5  ] • 
The  mechanism  is  basically  analogous  to  the  Gunn  effect  in  momentum  space. 
As  a  matter  of  fact,  when  the  Gunn  effect  was  discovered  in  1963.  the 
possibility  of  a  similar  effect  in  real  space  was  on  the  minds  of'  several 
scientista[s].  The  actual  realization,  however,  had  to  wait  for  the 
evolution  of  the  molecular  beam  epitaxy  (M3E)  technique  and  the  concept  of 
modulation  doping1]!  ,7] .  The  real  3pace  transfer  mechanism  is  based  on  the 
thermionic  emission  of  hot  electrons  from  the  high  mobility  GaAs  layer  to 
the  low  mobility  AlGaAs  layer  in  the  presence  of  a  high  electric  field. 
This  transfer  leads  to  a  negative  differential  resistance  (NDR)  .just  as  in 
the  Gunn  effect.  The  advantage  of  the  real  space  transfer  mechanism, 
however,  lies  in  its  degree  of  control  of  device  characteristics.  Negative 
differential  effects  are  important  for  their  possible  applications  to 
microwave,  switching,  and  memory  devices.  The  present  work  also  aims  to 
contribute  generally  to  the  theory  of  high  field  transport  in  semiconductor 
heterostructures.  Additionally,  the  theoretical  techniques  used  to  study 
this  effect  gives  us  the  understanding  and  the  tools  to  attack  another 

interesting  problem  of  semiconductor  physics  -  the  rigorous  treatment 

of  impact  ionization  in  semiconductors. 


Impact  ionization  is  a  very  important  hot  carrier  effect.  It  directly 
influences  the  performance  of  avalanche  phc todetectors]3 ]  and  IHPATT 
diodes]?].  It  also  determines  the  ultimate  performance  limits  of  small 


semiconductor  devices. 
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and  FETs'  111.  Nevertheless, 
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previous  theories  of  Impact  ionization  have  b*en  limited  by  various 
assumptions  and  also  by  the  use  of  several  adjustable  parameters,  such  as 
the  scattering  mean  free  path.  The  effective  mass  approximation  and  small 
non-parabolicity  correction  have  been  used  in  spite  of  the  fact  that 
carriers  gain  more  than  1  eV  in  energy  (as  measured  from  the  bottom  of  the 
conduction  band) .  In  our  work  we  have  attempted  to  include  a  realistic 
band  structure  in  the  Monte  Carlo  simulation  of  high  field  transport  in 
GaAs.  This  method  can  calculate  not  only  the  impact  ionization  rate  but 
also  other  quantities  of  interest,  such  as  drift  velocity  and  distribution 
functions.  The  inclusion  of  a  realistic  band  structure  is  expected  to 
expand  our  understanding  of  high  field  transport  to  extremely  high  electric 
fields  . 

In  Chapter  2  the  real  space  electron  transfer  mechanism  is  analyzed 
using  the  electron  temperature  model  and  the  Monte  Carlo  technique.  Then, 
various  effects  which  are  neglected  in  the  analysis  are  described.  Some 
comments  will  also  be  made  on  the  application  of  this  mechanism. 

In  Chapter  3  a  new  Monte  Carlo  method  using  a  realistic  band  structure 
is  described.  The  method  is  applied  to  the  study  of  high  field  transport 
in  GaAs.  On  this  basis,  a  detailed  study  of  the  band  structure  dependence 
of  impact  ionization  in  GaAs  is  given. 
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2.  REAL  SPACE  ELECTRON  TRANSFER  IN 
GaAs-Al  jjai  _xAs  HETSROSTRUC TUBES 

2 .1  General  Discussion 

The  basic  structure  consists  of  alternating  GaAs-Al  XA3  layers,  or 
other  appropriate  lattice-matched  materials  with  dissimilar  band-gap 
energies  and  carrier  mobilities.  Typical  layer  dimensions  are  from  100  8 
to  a  few  tenths  of  a  micron.  Such  structures  have  been  realized  by 
utilizing  molecular  beam  epitaxial  (MBE)[l  2-1  4],  vapor  phase  epitaxial 
(VPS ) [l 5 » 1 o] ,  or  liquid  phase  epitaxial  (LPE)[l7,13]  growth  techniques.  In 
this  work  we  are  concerned  with  GaAs-AlxGa^  As  multilayers.  For  0  <  x  < 
0.45  the  r direct  gap  of  A1  Ga^_xAs  changes  as[l9] 

Er  =  1.424  +  1. 247x  (eV) .  (2.1) 

§ 

In  GaA3-A.lxGa^  xAs  hetsrostruc tures  approximately  33  o  of  this  band  gap 
discontinuity  is  in-  the  conduction  band[20].  The  potential  barriers  form  a 
rectangular  potential  well  for  electrons.  At  thermal  equilibrium  electrons 
reside  at  the  minimum  of  the  potential  wells,  i.e.,  in  the  GaAs  layers. 
Application  of  a  high  electric  field  parallel  to  the  layer  interfaces  of 
thi3  structure  will  result  in  heating  of  the  electrons.  When  the  mean 
kinetic  energy  of  the  electrons  becomes  comparable  to  the  potential  barrier 
height,  AE,  they  can  be  thermionically  emitted  into  the  AlGaAs.  If  the 
mobility  of  electons  in  the  GaAs  layer  is  much  higher  than  the  mobility  in 
the  AlGaAs  layer,  the  sample  should  exhibit  negative  differential 
resistance  (NDR). 
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To  properly  describe  real  space  electron  transfer,  it  is  necessary  to 
account  for  the  scattering  mechanisms  encountered  by  electrons  being 
transported  at  high  kinetic  energies  in  the  potential  well.  A  mechanical 
analogue  to  real  space  transfer  is  provided  by  the  example  of  a  ball 
rolling  down  a  chute.  The  ball  will  3tay  in  the  chute  if  its  kinetic 
energy  remains  small.  However,  if  the  ball  gains  adequate  kinetic  energy, 
then  an  obstacle  can  scatter  the  ball  out  of  the  chute.  A  similar  effect 
occurs  in  a  layered  heterostructures,  where  electrons  drift  in  a  potential 
well  under  the  influence  of  a  high  electric  field. 

In  order  to  realize  large  difference  of  mobilities  in  two  materials 
modulation  doping[?,22]  of  the  layers  may  be  utilized.  Fig.  2.1 
schematically  shows  the  band  structure,  doping  distributions,  and  electron 
mobility  of  modulation-doped  GaAs-AlxGa1  xAs  layers.  The  AlGaAs  layers  are 
intentionally  doped  to  a  density  of  10  -10  cm~  ,  whereas  the  GaAs 

layers  contain  only  unintentional  background  impurities  (lO1^  ~  lO1^ 
cm-^).  At  thermal  equilibrium  the  electrons  reside  in  the  GaAs  layers. 
When  separated  more  than  200  R  these  electrons  experience  strongly  reduced 
impurity  scattering.  Hess  has  shown  that  the  differential  scattering  rate 
due  to  ionized  impurities  decreases  exponentially  with  spatial 
distance[23].  This  mobility  enhancement  by  spatial  separation  of  electrons 
from  donors  wa3  first  proposed  by  Ssaki  and  Tsu[4],  and  later  verified 
experimentally  by  Dingle  et  al.[7].  They  have  observed  a  mobility  of  5000 
cm^/Vsec  at  300  K,  and  even  higher  mobilities  (~  15000  cm^/Vsec''  at  low 
temperatures.  The  AlGaAs  layers  can  be  made  strongly  compensated,  and 
thereby  the  mobility  in  these  layers  can  be  very  low  (_  500  cm^/Vsec  or 
less)  .  In  this  manner  a  mobility  ratio  of  ten  or  more  C3n  be  realized  in 
modulation  dt^ed  structures. 


Energy  of 
Conduction 

Band  Edge  Ec  (meV)  Mobility  /z  (cm2A 


Fig.  2.1.  Schematic  diagrams  of  doping  density,  electron  mobility,  and 
conduction  band  energy  of  modulation-doped  layers. 
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Sines  the  carrier  temperature  i3  determined  by  the  power  balance 
equation: 


epF 


2 


dE 

dt 


(2.2) 


and  the  energy  loss  rates  are  similar  in  the  two  materials,  the  heating  of 
carriers  occurs  only  in  the  GaAs  layer.  These  hot  electrons  transfer  in 
space  to  the  adjacent  AlGaAs  layers  and  then  thermal ize.  Some  of  the 
cooler  electrons  in  the  AlGaAs  layers  can  undergo  reverse  transfer  into  the 
GaAs  layer.  In  order  to  take  account  of  this  energy  and  momentum  exchange 
between  the  GaAs  and  the  AlGaAs  layers,  one  must  solve  the  3oltzmann 
equation: 


3_f 
3 1 


->■  ->■ 
V,  f  +  v 
k 


(2.3) 


with  appropriate  boundary  conditions.  In  addition  to  the  spatial 
inhomogeneity  of  this  system,  there  are  several  other  factors  which 
complicate  the  analysis: 

i)  At  high  electric  fields  the  energy  distribution  functions  in 
polar  semiconductors  are  highly  non-Maxwellian. 

ii)  There  is  an  energy  exchange  between  the  layers  because 
energetic  carriers  are  flowing  out  of  the  GaAs  and  cold 
carriers  are  returning. 

iii)  Electron-electron  interactions  must  be  taken  into  account, 

iv)  The  effects  of  higher  L  minima  in  two  marerials  should  be 
considered . 

v)  There  are  several  other  effects  which  make  the  analysis 


extremely  difficult,  3uch  as  two  dimensional  effects,  band 
bending,  reflection  and  trapping  of  electrons  at  the 
interfaces,  etc. 

To  thi3  end  we  have  chosen  two  methods;  the  electron  temperature  model  and 
the  Monte  Carlo  method.  The  electron  temperature  model  assumes  nearly 
Maxwellian  distributions  in  the  two  materials.  The  Monte  Carlo  method  is  a 
one  particle  simulation  of  electron  transport  in  a  simplified  potential 

profile.  Each  method  applies  only  under  certain  conditions.  Detailed 
discussion  of  the  two  models,  assumptions,  and  the  results  of  the 

calculations  are  presented  in  the  next  two  sections. 

In  both  methods  we  have  deliberately  chosen  typical  layer  widths 
larger  than  400  S.  This  avoids  the  complications  arising  from  size 
quantization  effect  and  other  two  dimensional  effects.  Some  of  the  two 

dimensional  effects  are  discussed  in  Section  2.4.  The  choice  of  large 

layer  widths  also  allows  the  use  of  the  semiclassical  Boltzmann  equation. 
Otherwise  one  must  resort  to  a  quantum  mechanical  treatment.  With  large 
layer  widths,  the  potential  fluctuations  due  to  the  statistical 
distribution  of  the  impurities  can  also  be  neglected.  Other  simplifying 
assumptions  in  the  models  are  taken  up  and  reexamined  in  Section  2.5. 
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2 .2  Slactron  Temperature  Model[24 ] 

2.2.1  Description  of  the  model 

In  polar  semiconductors  the  energy  distribution  of  electrons  can  be 
highly  non-Maxwellian,  since  polar  optical  scattering  is  inelastic  and  i3 
not  randomizing  (see  Fig.  A. 3  in  Appendix  1).  To  make  the  analysis 
explicitly  tractable,  however,  we  assume  a  Maxwellian  form  for  the 
isotropic  part  of  the  distribution  function.  Thi3  assumption  is  partly 
justified  by  the  fact  that  the  electron  density  is  high  in  the  GaAs  layer 
at  the  3tart  of  the  electron  transfer,  3ince  the  0a As  layer  collects 
electrons  from  the  neighboring  AlGaAs  layers.  Typically  the  GaA3  layer 
will  have  an  electron  density  of  1018  „  10  1  ^  cm"^ .  Electron-electron 
collisions  will,  therefore,  randomize  the  energy  gained  in  the  electric 
field  direction  and  establish  a  Maxwellian  distribution.  Above  the  band 
edge  of  the  AlGaAs  the  electrons  in  the  GaAs  will  follow  a  "cooler" 
distribution  because  of  the  reverse  transfer  of  cooler  electrons  from  the 
AlGaAs.  This  will  be  shown  as  a  result  of  Monte  Carlo  simulations  in 
Section  2.3*3  (Fig.  2.25).  However,  since  the  GaAs  layer  is  sufficiently 
narrow,  we  can  use  a  single  electron  temperature  in  the  GaAs. 

If  there  were  no  interaction  between  the  GaAs  and  AlGaAs  layers,  the 
electron  temperature  in  the  AlGaAs  layer  would  be  constant  and  close  to  the 
lattice  temperature,  T  .  However,  this  is  not  a  realistic  assumption.  Due 
to  the  transfer  of  hot  electrons  from  the  GaAs  layer,  electrons  close  to 
the  GaAs  layer  will  have  considerably  higher  temperature  than  the  lattice 
temperature.  As  electrons  move  away  from  the  interface  they  thermalize  to 
the  equilibrium.  Therefore,  we  need  to  consider  the  position  dependence  of 
the  electron  temperature  in  order  to  account  for  the  power  flow 
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perpendicular  to  the  interface. 

Fig.  2.2  shows  the  model  used  in  the  calculations.  We  assume  a 
position-dependent  electron  temperature,  T^(x) ,  and  a  position-dependent 
quasi-Ferni  level,  E^Cx),  in  the  AlGaAs  layer.  Similar  methods  have 
recently  been  employed  by  other  workers  to  analyze  thermionic  emission  in 
metal-GaAs  contactsf 25]*  Both  the  concept  of  the  electron  temperature  and 
of  quasi-Fermi  levels  require  high  carrier  densities  in  order  to  be  valid. 
In  the  thin  GaAs  layer  we  assume  a  position- independent  electron 
temperature,  Tq>  and  quasi-Fermi  level,  E^  .  With  this  assumption  we  only 
need  to  solve  the  position  dependence  of  the  electron  temperature  and  the 
quasi-Fermi  level  in  the  AlGaA3  layers.  The  thermal  conduction  of  hot 
electrons  from  the  GaAs  layer  into  the  AlGaAs  layers  is  then  accounted  for 
by  the  slope  of  T^(x)  at  the  interface  boundary.  Although  at  the  boundary 
the  slope- of  T  in  the  GaAs  layer  is  zero  (since  we  have  assumed  T  to  be 
constant)  ,  we  assume  that  the  same  amount  of  energy  as  given  by  the  slope 
of  Tt(x)  is  extracted  from  the  GaAs. 

jJ 

2.2.2  Calculations 

We  proceed  to  assume  a  distribution  function  of  a  form: 

f  =  f  (  e,  x)  +  g(  e,  x)k  +  h(  e,  x)k 

O  J  wv 

where  f  ( z,x)  is  a  Maxwellian  distribution: 

[Ep2  ^  -  £ 

f o  (  e  ’  x)  =  exP|_  kTL(x)" 
in  the  AlGaAs  layer.  The  coordinate  system  is  shown,  in  Fig.  2.2.  The 


(2.4) 


(2.5) 


electric  field  is  in  the  y-direction.  The  second  tern  of  Eq.  (2.4) 
represent  the  drift  term  of  the  distribution  along  the  electric  field  as  is 
usually  used  for  a  spatially  uniform  system'] 26].  The  third  term  is  to  take 
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into  account  the  spatial  nonunifo rmity.  It  represents  the  displacement  of 
the  distribution  perpendicular  to  the  electric  field.  Thi3  term  is 
essential  to  account  for  the  energy  flow  normal  to  the  field  direction. 
The  isotropic  part  of  the  distribution  function  in  the  GaAs  layer  is  given 
by  a  form  similar  to  Sq.  (2.5)  in  terms  of  Tg  and  E?1  ,  but  without  the 
position  dependence. 


As  is  normally  done,  we  then  separate  the  variables  into  equations  for 
the  spatially  symmetric  part  and  the  drift  terms[26].  The  method  of 
moments  [27]  is  then  utilized  to  obtain  the  differential  equations  for  Tt(x) 
and  379(x).  The  set  of  equations  obtained  is  given  below: 


2ho)  ^/2  exp (Z  -Z  )-l  .  .. 

- -|  eE  - %  y-j -  ZTi/Z  exp ( Z  / 2 )  K  (Z. /2) 

•  1  o  exp(Z  )-l  L  L  oL 


mm. 


=  ey2F  +  2e 


ry  _  O  /  y  J  ^ 

—  (kTT)  exp(-EF2/kTL) 


dx 


[(kTL> 


7/2 


exp(EF7/kTL) 


(2.6) 


and 


2  - 

^2  i"(kTL)5/2  exp(EF2/kTL)]  =  0 
dx 


(2.7) 


kTL  kT 

where  ZT  =  r —  ,  Z  =  r— ' -  T  i3 

L  ftw  o  ho;  o 

o  o 

mobility  in  the  AlGaA3  layer,  7  is 
the  Bessel  function  of  the  zeroth 
3T7(x)  and  T^(x)  have  been  omitted 
Eq.  (2.6)  is  the  rate  of  energy  loss 


the  lattice  temperature,  U2  is  the 
the  applied  electric  field,  and  X  is 
order.  The  position  dependences  of 
in  the  notations.  The  left  side  of 
due  to  polar  optical  scattering1"  26^ , 
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and  the  first  term  on  the  right  side  is  the  power  input  from  the  applied 
electric  field.  The  second  tern  on  the  right  represents  the  energy  flow 
from  the  GaAs  due  to  the  transfer  of  hot  electrons.  Applying  energy 
conservation  at  the  boundary,  we  then  obtain  a  similar  power  balance 
equation  relating  Tg  and  F  in  the  GaAs  layer; 


eE 

o 


1/2 

exp(Zo)-l  e 


exp (Z  /2)  K  (Z  /2) 
e  o  e 


5u,  ,  dT  (0) 

'  *“1F  +  k  V0)  -3T- 


(2.8) 


kT 

where  z  =  r — and  u,  is  the  mobility  in  the  GaAs  when  only  polar  optical 
e  noi  1  i 

o  ! 

scattering  is  operative.  This  mobility  is  given  by; 


u 


1 


exp(-Ze/2) £[  exp(Zo-Ze)+l  ]  K^/2) 


+  [  exp (Z  -Z  )-l  ]  K  (Z  / 2 )T- ~ 1 
r  o  e  o  e  J 


(2.9) 


where  is  the  Bessel  function  of  the  first  order,  and  N  is  given  by 


N  = 


q  exp(hcu  /k  T  )-l 
o  b  ° 


(2.10) 


The  second  term  on  the  right  side  of  Eq.  (2.3)  represents  the  power  flowing 
out  of  the  GaAs  into  the  AlGaA3  layers.  1^(0)  and  dTL(0)/ax  are  evaluated 
at  the  boundary  x  =  0.  We  have  used  the  fact  that  this  power  flow  is  due 
only  to  those  electrons  in  the  GaAs  with  energy  higher  than  the  AlGaAs  band 
edge.  Combining  Eqs.  (2.6)  and  (2.7),  we  can  eliminate  E^0(x)  and  obtain 
the  differential  equation  in  terms  of  T  (x)  only: 

±J 
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d'TL(x) 

dx2 


2e 


5k  T^(x) 


0  2hw 
F2- 


exp(l-zL) 


-1 


,1/2 


/  *  o  exp  ( Z  )  —  1  L 

M2'.'1Tm2  L 


x  exp  (Z  / 2 )  K  (Zt/2) 
L  o  L 


(2.11) 


Thi3  equation  is  now  solved  numerically  with  the  boundary  condition;  Tt(q) 
■  Te,  and  dTL(x)/dx  =  0  at  TL  =  TQ. 


To  determine  the  quasi-Femi  levels,  E^.  and  Et,~(x),  we  need  two 
additional  conditions.  One  is  obtained  from  the  condition  j.  =  0.  ‘Ve 
simplify  this  condition  to  the  balance  of  the  thermionic  currents  in  both 
directions  at  the  interface.  Under  collision- free  transport  conditions  the 
thermionic  current  from  the  3a  As  to  the  AllaAs,  is  given  after 


3ethe[ 23]  by; 


* 


em 


'1-2 


2tt 


1  (kT  )2 
e 


exp 


-ef1-ae- 


kT 


(2.12) 


The  current  flowing  from  the  AlGaAs  to  the  Oa As  layer  is; 


em 


‘2-1 


„  rE  2-AE 
(kT  )  “  exp! — rp 


2tt 


2h3  L  kT 


(2.13) 


where  T^  is  evaluated  at  x  =  0.  Under  3teady-state  conditions  we  have  .j,  0 
=  Jg-i  *  The  other  additional  condition  arises  from  conservation  of  the 
total  number  of  electrons,  N  =  N 1  +  y  .  and  are  calculated  using  the 

electron  temperature  and  quasi-Ferai  level  in  each  layer. 
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2.2.3  Results 

Prom  3q.  (2.3)  we  obtain  the  variation  of  electron  temperature,  T 

*■  e 

with  the  electric  field  in  the  GaAs  layer.  The  result  is  3hown  in  Fig.  2.3 
with  the  mobility  in  the  AlGaAs  layer  as  a  parameter.  The  dashed  curve  is 
for  bulk  GaAs,  with  only  polar  optical  scattering  operative.  This 
corresponds  to  treating  the  GaA3  and  the  AlGaAs  layers  independently  (in 
other  words,  =  0).  As  is  well  known[29> 30 ] ,  at  an  electric  field  of 


about  3-4  kV/cm 

the 

electron 

temperature 

goes 

to  infinity 

and  "polar 

runaway"  occurs. 

When  the  t 

ransport  of 

hot 

electrons  is 

;aken  into 

account,  however, 

the 

electron 

temperature 

does 

not  increase 

as  rapidly. 

The  rate  of  increase  is  slower  for  higher  mobility  in  the  AlGaAs,  because 
of  the  larger  power  out-flow.  This  situation  is  quite  analogous  to  the 
Gunn  effect[30, 31  ]  where  intervalley  scattering  supplies  an  additional 
energy  loss  mechanism[32] . 

Fig.  2.4  shows  the  position  dependence  of  both  the  electron 
temperature,  T^(x)  ,  and  the  Fermi  level,  E-^Cx)  ,  for  everal  values  of 
electric  field.  The  dashed  line  indicates  the  GaAs  potential  well  of  250 
meV  depth.  The  interface  i3  set  at  x  =  0.  It  can  be  seen  that  both  T  (x) 
and  Ej^x)  approach  their  equilibrium  values  within  a  distance  on  the  order 
of  0.1  p.  The  small  spikes  in  the  quasi-Fermi  levels  which  occur  at  the 
boundary  arise  from  the  Imposed  condition  of  thermionic  current  balance  and 
from  the  small  difference  in  effective  masses.  These  spikes  are  of  no 
physical  significance.  The  rapid  decrease  of  the  quasi-Fermi  level  in  the 
GaAs  with  higher  fields  indicates  the  depletion  of  electrons  in  this  layer, 
and  hence  the  transfer  of  electrons  into  the  AlGaAs.  The  variation  of  the 
fraction  of  electrons  in  the  GaAs  layer  -with  increased  field  is  illustrated 


Electron  Temperature,!", 
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Fig.  2.3.  Variation  of  electron  temperature  in  the  GaAs  layer  with 

applied  field  for  various  mobilities  in  the  A1  Ga,  As  laver. 

x  1-x 
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in  Pig.  2.5*  The  transfer  is  larger  for  smaller  mobility  values,  a 
consequence  of  the  greater  carrier  heating  for  smaller  mobility  a3  seen  in 
Fig.  2.3*  Here  again  the  analogy  with  the  Sunn  effect  should  be  noted[32]. 

The  current-voltage  characteristics  are  straight  forward  to  calculate 
from  this  model.  The  results  of  the  calculations  are  shown  in  Fig.  2.6. 
The  magnitude  of  the  NDR  changes  dramatically  with  the  mobility  in  the 
AlGaA3  layer.  For  purposes  of  comparison  we  include  in  Fig.  2.6  (insert) 
the  velocity-field  characteristics  calculated  by  Fawcett  et  al.[33]  in  1970 
for  the  Sunn  effect.  The  parameter  for  their  curves  was  the  intervalley 
deformation  potential.  Since  no  reliable  data  on  the  strength  of 
intervalley  scattering  wa3  available  in  1970,  the  intervalley  deformation 
potential  was  used  as  a  parameter. 

2.2.4  Comparisons  -with  the  Sunn  effect 

The  two  sets  of  curves  in  Fig.  2.6  show  surprisingly  similar  features. 
Actually  our  discussion  of  real  space  transfer  almost  parallels  that  of  the 
Gunn  effect,  but  with  two  crucial  differences.  First,  in  our  mechanism 
electrons  leave  the  high  mobility  GaAs  layer  by  thermionic  emission  and  are 
transferred  in  real  space  to  the  low  mobility  .41 GaAs  layer.  In  the  Gunn 
effect,  on  the  other  hand,  electrons  transfer  from  one  valley  of  high 
mobility  to  another  of  low  mobility  in  momentum  or  k-space.  Second  and 
most  Important  is  the  fact  that  our  device  characteristics  can  be 


controlled  to  a  greater 

degree 

than  with  a  device 

utilizing 

^  Vi a 

Gunn 

effect.  As  already  seen 

in  Fig , 

.  2.5  (insert),  the 

curves  for 

f  Vjo 

Gun  n 

effect  have  as  a  parameter 

the  int 

ervalley  deformation 

po  tential , 

wh  lc 

h  is 

a  material  property  and  cannot  be  changed.  The  intervalley  deformation 
potential  between  the  r  and  L  valleys  i3  now  known  to  be  around  1  x  IT 


Fraction  of  Electrons  in  the  GaAs 


Current,  J  (arbitrary  units ) 


0  5  10  15 

Electric  Field,  F  (kV/cm)  L?.l541 


Fig.  2.6.  Current-voltage  characteristics  for  a  real-space  transfer  device 
with  various  values  of  electron  mobility  in  the  A1  Ga.  As 
layer.  The  inset  shows  the  velocity-field  characteristics  of 
the  Gunn  effect  in  GaAs  for  various  intervalley  deformation 
potentials  [33  ]  . 
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eV/cm  in  GaAs[34]«  On  the  other  hand,  the  parameter  in  Fig.  2.6  for  the 
real  3pace  transfer  process  is  the  mobility  of  the  AlOaAs  layer,  which  can 
be  controlled  by  adjusting  the  doping  of  the  layer.  The  characteristics 
can  also  be  changed  by  varying  other  device  parameters,  3uch  as  the  layer 
dimensions  and  the  potential  barrier  height. 

2.2.5  Influence  of  material  parameters  on  the  real  space  transfer 

The  potential  barrier  height  can  be  controlled  by  changing  the  A1  mole 
fraction  of  the  AlOaAs .  The  effect  of  this  parameter  on  the 
current-voltage  characteristics  is  shown  in  Fig.  2.7.  For  this  particular 
set  of  parameters  the  threshold  field  for  the  onset  of  NDR  can  be  varied 
between  2  and  3  kV/cm  and  the  peak-to- valley  ratio  between  1.3  and  2.0. 
Fig.  2.3  shows  the  fraction  of  electrons  in  the  OaAs  for  the  same 
parameter.  The  transfer  of  electrons  is  more  abrupt  for  the  deeper 
potential  well,  which  results  in  larger  peak-to- valley  ratio  as  shown  in 
Fig.  2.7.  This  occurs  because  the  deeper  the  potential  well  is,  the  more 
the  electrons  can  be  heated  before  they  begin  to  transfer,  and  hence  the 
larger  gradient  of  carrier  temperature  results  at  the  interface  at  the 
beginning  of  transfer.  In  this  electron  temperature  model  we  do  not  take 
into  account  the  effects  of  the  conduction  band  L  minima,  which  are  known 
to  be  located  approximately  330  meV  above  the  T  minimum  of  GaAs[34, 36] .  We 
have  deliberately  chosen  values  of  AE  such  that  the  L  minima  are  above  the 
band  edge  of  the  AlOaAs.  The  effects  of  the  L  minima  will  be  discussed 
more  extensively  in  Section  2.3- 

The  layer  dimensions  can  be  varied  by  changing  the  crystal  growth 
parameters.  Fig.  2.9  shows  the  electron  temperature  in  the  OaAs  as  a 
function  of  electric  field  for  several  thicknesses  of  the  OaAs  layers. 
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Fig.  2.7.  Current-voltage  characteristics  for  various  potential  barrier 
heights.  Note  that  the  GaAs  L  mimina  are  above  the  conduction 
band  edge  of  the  A1  Gan  As  in  the  cases  shown. 


Electric  Field,  F  (kV/cm) 


10 


L' 


Fraction  of  electrons  in  the  GaAs  layer  as  a  function  of 
electric  field  for  various  potential  barrier  heights. 
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2.9.  Variation  of  electron  temperature  in  the  GaAs  layer  with 
applied  field  for  various  thicknesses  of  the  GaAs  layer. 
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degree  of  carrier  heating  is  smaller  for  thinner  GaAs  layers  because  of  the 
larger  cooling  effect  of  the  surrounding  AlGaA3  layers.  Fig.  2.10  shows 
the  effect  of  changing  the  GaAs  thickness  on  the  current-voltage 
characteristics. 

The  lattice  temperature  dependence  of  the  current- voltage 
characteristics  is  another  interesting  subject[35].  Fig.  2.11  shows  the 
result  of  the  calculation  in  the  temperature  range  77  ~  400  K.  In  this 
calculation  the  electron  mobility  in  the  AlGaAs  layer  is  assumed  to  be 
constant  at  100  cm'Vvsec.  In  heavily  doped  materials,  and  in  particular, 
compensated  materials,  the  mobility  is  known  to  vary  very  little  with 
tempera ture[ 37] .  Since  the  mobility  in  the  GaAs  layer  increases 
considerably  a3  the  temperature  decreases,  the  larger  psak-to-valley  ratio 
results  at  lower  temperature.  According  to  our  calculation,  the 
peak- to- valley  ratio  improves  from  1.56  at  400  X  to  7.75  at  77  X.  This 
suggests  that  the  experimental  manifestation  of  the  real  space  transfer 
mechanism  can  be  more  easily  realized  at  lower  temperature. 

Although  the  electron  temperature  model  is  based  on  several 
simplifying  assumption,  some  of  which  are  difficult  to  justify,  it  provides 
a  qualitative  method  with  which  to  analyze  the  real  3pace  transfer 
mechanism.  The  model  is  exact  in  the  high  carrier  density  limit.  The 
Monte  Carlo  method  is  easily  tractable  in  the  other  extreme,  i.e.  the  low 
density  limit.  As  will  be  seen  the  results  obtained  from  both  method  show 
surprisingly  similar  features  (such  as  dependence  on  barrier  height,  AlGaAs 
mobility,  etc.;.  Thi3  suggests  that  these  features  are  in  fact  general  and 
do  not  depend  on  the  assumptions  of  the  model. 
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- -3  Monte  Carlo  3 imulation  of  Real  Space  Electron  T ransfer 
2.3*1  Introduction 

The  Monte  Carlo  method[3S]  used  in  transport  theory  i3  based  on  a 
computer  simulation  of  the  carrier  motion  of  an  individual  electron  in  an 
electric  field.  Physical  quantities  of  interest  (such  as  drift  velocity) 
are  obtained  as  proper  statistical  averages  among  sample  members  of  an 
ensemble.  The  Monte  Carlo  procedure  has  been  shown  to  yield  a  distribution 
function  which  is  a  solution  of  the  3oltzmann  equation[33]  and  thus 
provides  a  powerful  method  to  solve  the  3oltzmann  equation  under  quite 
general  conditions.  It  has  been  applied  with  great  success  to  the  scudy  of 
high  field  transport,  and  in  particular,  to  the  Cunn  effect[33]. 

The  Monte  Carlo  method  follows  a  carrier  experiencing  successive  drift 
and  scattering  event.  The  scattering  event  and  the  duration  of  each  drift 
are  determined  by  random  numbers.  The  simulation  is  repeated  a  sufficient 
number  of  times  to  minimize  statistical  fluctuations. 

The  effectiveness  of  the  Monte  Carlo  method  lies  in  its  simple 
principle  and  its  flexibility.  The  inclusion  of  various  complex  scattering 
mechanisms,  time  dependence,  and  physical  boundaries  in  position  space  can 
be  done  with  minor  elaboration  of  the  algorithm.  It  is  this  flexibility 
that  makes  the  Monte  Carlo  method  especially  useful  for  the  simulation  of 
real  space  transfer. 
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2.3*2  Simulation  procedure 

Fig.  2.12  shows  the  flow  chart  for  the  Monte  Carlo  simulation  of  real 
space  transfer  in  GaAs-AlGaAs  heterostructures.  The  quantities  r,  £,  and  E 
denote  the  position  vector,  the  k-vector,  and  the  energy  of  the  electron 
respectively.  The  appearance  of  r^  in  the  flow  chart  means  that  a  random 
number  'uniformly  distributed  between  0  and  1  is  necessary  at  that  stage. 

The  lower  right  portion  of  the  flow  chart  is  specific  to  the 
simulation  of  real  space  transfer  mechanism.  The  spatial  configuration  is 
shown  in  Fig.  2.13.  The  model  assumes  an  abrupt  potential  barrier  of 
height  A3  in  the  GaAs.  3and  bending  effects  are  neglected  as  in  the 
electron  temperature  model.  The  ba3ic  cell  from  x  =  0  to  x  =  d^  +  i  i3 
cyclically  repeated  in  the  simulation.  The  calculation  of  the  transmission 
coefficient  is  described  in  detail  in  Section  2.5*2. 

All  the  formula  and  material  constants  necessary  for  the  calculation 
of  the  scattering  rates  are  given  in  Appendices  1  and  2.  After  the 
scattering  rate  of  each  mechanism  has  been  calculated,  the  total  scattering 
rate  i3  obtained  by  summing  over  all  the  possible  scattering  mechanisms  at 
each  value  of  electron  energy.  The  probability  table  is  then  constructed, 
which  is  used  in  the  determination  of  the  scattering  machanism.  This  table 
contains  the  relative  probability  of  each  scattering  mechanism  at  every 
value  of  electron  energy.  The  total  scattering  rate  is  stored  in  memory  in 
the  form  of  the  3lope  and  the  intercept  which  linearly  interpolate  the  rate 
oetween  two  neighboring  energy  points.  This  table  of  the  slope  and  the 
intercept  is  used  in  the  look-up  scheme  to  calculate  the  total  scattering 
rate  at  an  arbitrary  electron  energy[39]>  The  integration 
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The  simulation  begins  by  releasing  an  electron  with  thermal  energy 
from  the  center  of  the  GaAs  well  and  in  a  randomly  selected  direction.  The 
electron  then  undergoes  scattering  interactions.  The  scattering  mechanism 
considered  in  the  calculations[54]  include  acoustic  phonon  scattering, 
optical  phonon  scattering,  piezoelectric  scattering,  equivalent  and 
non-equivalent  intervalley  scattering,  ionized  impurity  scattering,  and 
random  potential  alloy  scattering[73 1  -in  the  AlGaAs.  In  some  of  the 
simulations  only  polar  optical  scattering  and  intervalley  scattering  in  the 
GaAs  and  ionized  impurity  scattering  and  polar  optical  scattering  in  the 
AlGaAs  are  used  to  minimize  the  computer  time.  This  does  not  affect  the 
essential  features  of  the  results,  since  these  scattering  mechanisms  are 
the  most  dominant  ones.  Some  comparisons  will  be  made  between  these  two 
cases  in  the  next  section.  In  all  of  the  simulations  only  the  T  (000)  and 
Li(  1  11  )  conduction  bands  are  considered  for  GaA3,  and  only  the  F  (000 ) 
conduction  band  is  considered  in  AlGaAs. 

A  typical  simulation  consists  of  approximately  100,000  to  250,000 

interactions  to  obtain  accurate  velocity  estimates.  The  average  irift 

velocities  in  the  individual  materials  and  the  average  velocity  in  the 

<*  ~\ 

heterostructures  are  calculated  using  standard  velocity  estimators'. S3  J . 
For  the  ohmic  range  (in  the  GaAs  at  low  field  and  in  the  low-mobility 
AlGaAs)  the  diffusion  coefficient  is  calculated  first  from  a  maximum 


likelihood  estimate: 
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<  (y  -  vdt)2  >  =  2Dt  (2.15) 

where  v^  i3  the  drift  velocity,  and  t  is  the  duration  of  sampling.  Then 
the  Rinstein  relation  is  used  to  calculate  the  mobility.  This  method  has 
been  found  to  reduce  the  uncertainty  in  the  calculated  drift  velocity  at 
low  electric  field[40, 41 ] • 

2.3*3  Results 

Rig.  2.14  shows  an  example  of  the  electron  trajectory  in  real  space 
for  low  (2  kV/cm)  and  high  (6  kV/cm)  electric  fields  obtained  by  the  Monte 
Carlo  method.  The  narrow  region  between  two  lines  represents  the  CaAs 
layer,  whereas  the  region  outside  are  the  AlGaAs  layers.  At  low  electric 
field  (2  kV/cm)  the  electron  stays  inside  the  potential  well  as  it  drifts 
along  the  well,  since  the  average  electron  energy  ( „  0.06  eV)  is  much 
smaller  than  the  potential  barrier  height  (0.2  eV).  The  mean  free  path  of 
scattering  for  this  electric  field  in  undoped  GaAs  is  approximately  600  X. 
Therefore  the  electron  i3  reflected  by  the  interface  more  often  than  it  is 
scattered  by  a  phonon.  In  our  simulation  we  have  assumed  a  specular 
reflection  at  the  wall  with  no  loss  of  velocity  in  the  field  direction.  In 
reality  the  interface  scattering  may  play  an  important  role  in  determining 
the  mobility  in  the  GaAs  potential  well.  When  the  electric  field  is  high 
(6  kV/cm)  the  electron  gains  enough  kinetic  energy  ( ~  0.1R  eV)  so  that  it 
can  be  scattered  out  to  the  AlGaAs  layer  as  shown  on  the  right  side  of 
Fig.  2.14.  Since  the  electron  mobility  in  the  AlGaAs  is  much  lower,  the 
electron  moves  much  slower  in  thi3  layer,  and  the  mean  free  path  is  much 
shorter.  There  i3  a  possibility,  of  course,  that  the  electron  will  move 


t  tii  =400 A  t 


AIGaAs  0 
d2  =  4000A 


Below  Threshold  Above  Threshold 
F=  2  kV/cm  F=  6  kV/cm 


-S-L6I5 


Fig.  2.14.  Trajectory  of  an  electron  in  real  space  for  electric  fields 
below  and  above  the  threshold  for  the  real  space  transfer 
mechanism. 


back  into  the  GaAs  potential  well  after  it  moves  in  the  AlGaA3  layer  for  a 
while.  This  is  also  illustrated  in  the  figure.  When  observed  over  a 
sufficiently  long  time  this  individual  electron's  motion  should  be 
representative  of  all  the  electrons  in  the  GaAs-AlGaAs  structure.  This 
provides  a  one  particle  picture  of  real  space  transfer. 

The  velocity-field  characterictic  of  a  real  space  transfer  (RST) 

device  obtained  by  the  Monte  Carlo  method  is  shown  in  Fig.  2.15*  This 

simulation  has  been  done  by  Glisson  et  al.[S4]  including  all  the  scattering 

mechanisms.  For  comparison  the  bulk  characteristics  of  the  GaAs  and  the 

\llaA3  ure  shown.  The  GaAs  has  no  ionized  impurities,  and  shows  the  usual 

»tnn  effect  at  1  kV/cm.  The  low-field  mobility  is  3000  cm^/Vsec.  The 
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•-•a As  : 3  3trongly  compensated  -with  an  impurity  density  of  1  x  10  cm 
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in:  a  free  electron  density  of  1  x  10  cm  .  This  gives  an  electron 
0  -  of  approximately  500  cm  /Vsec  in  thi3  material.  The  transport  in 

*  -  remains  essentially  ohmic  for  the  range  of  electric  field  shown 

.0  *  he  figure.  The  impurity  density  should  be  considered  as  a  parameter 
fcr  varying  the  mobility  in  the  AlGaA3  layer.  In  reality  it  may  not  be 
necessary  :o  dope  the  layer  with  this  high  impurity  density  if  there  is 
another  alternative  way  to  grow  a  low-mobility  material.  The  A1  content  in 
the  AllaAs  is  chosen  to  give  a  potential  barrier  of  0.2  eV  (x  =  0.18)  in 
the  3ST  device.  The  real  space  transfer  structure  shows  NDR  at  2.4  kV/cm 
with  a  peak  velocity  of  1.5  x  10^  cm /sec. 

For  comparison  we  show  in  Fig.  2.1o  the  calculated  velocity-field 
characteristics  -with  only  polar  optical  and  intervalley  scatterings  in  the 
GaAs  and  polar  optical  and  ionized  impurity  scatterings  in  the  AlGaAs.  In 
this  case  we  have  fewer  mechanisms  for  energy  and  momentum  loss.  As  3 
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result,  the  Gunn  effect  occurs  at  a  slightly  lower  electric  field.  The  low 
field  mobility  for  GaAs  and  the  RST  device  is  slightly  higher  for  the  same 
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reason.  The  usual  Brooks-Herring  formula  for  ionized  Impurity 
3cattering[27]  has  been  used  in  this  calculation,  whereas  Glisson  et  al . 
have  assumed  a  cut-off  distance  for  the  screened  Coulomb  interaction.  The 
Brooks-Herring  formula  gives  a  mobility  of  only  270  cm“/Vsec  for  a  lower 
impurity  density  of  4  x  10  °  cm-7  (and  the  same  free  electron  density  of  1 

4  <-r  -r 

x  10  '  cm-7).  The  reason  for  the  discrepancy  is  that  for  very  high  dopings 
the  average  distance  between  impurities  becomes  smaller  than  the  Debye 
screening  distance.  As  a  result,  the  Brooks-Herring  form -ala  (which  only 
considers  an  isolated  impurity)  overestimates  the  interaction.  For  the 
analysis  of  real  space  transfer,  however,  thi3  discrepancy  is  irrelevant 
because  the  characteristics  are  determined  mainly  by  the  mobility.  As  long 
as  the  mobilities  are  similar,  the  two  results  are  expected  to  be  similar. 

Fig.  2.17  shows  the  relative  numbers  of  electrons  in  the  GaAs  and  the 
AlGaAs  layers  using  the  same  parameter  values  as  shown  in  Fig.  2.15.  It 
can  be  3een  that  even  for  fields  well  above  threshold  a  substantial 
fraction  of  electrons  remains  in  the  GaA3.  These  residual  electrons  in  the 
GaA3  layer  might  hinder  the  application  of  this  mechanism  to  switching 
devices  3ince  they  constitute  a  considerably  large  current  in  the  "off" 
state.  The  ratio  of  electrons  in  the  two  materials  can  be  Improved, 
however,  by  reducing  further  the  mobility  in  the  AlGaA3.  In  Fig.  2.13  the 
mobility  is  reduced  to  about  50  cm^/Vsec  by  increasing  the  impurity  density 
to  1  x  10  cm-7.  In  thi3  case  less  than  10  %  of  the  electrons  remain  in 
the  GaA3  at  3  kv/cm.  The  peak  velocity  of  1.6  x  13^  cm/sec,  threshold 
field  of  2.3  kV/cm,  peak-to- valley  ratio  of  1,  and  a  negative  mobility 
magnitude  above  threshold  of  greater  than  1  x  ‘  0^  cm-,/7sec  are  very 
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attractive  parameters  for  Gunn-type  device  considerations.  On  the  other 
hand,  if  the  AlGaAs  mobility  is  increased,  a  substantial  number  of 
electrons  remain  in  the  GaAs.  Now  the  transport  property  of  the  structure 
is  determined  by  both  the  GaAs  and  the  AlGaAs.  This  is  illustrated  in 
Fig.  2.19  for  an  AlGaAs  mobility  of  4000  cm  /Vsec.  The  curve  shows  the 
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saturating  characteristic  and  a  peak  velocity  of  2  x  10  cm/sec  with  very 
small  negative  resistance.  This  illustrates  the  possibility  of 
artificially  creating  the  velocity- field  characteristics  using  a  layered 
heterostructure . 

As  in  the  electron  temperature  model  it  is  possible  to  study  the 
effects  of  changing  various  device  parameters  with  the  Monte  Carlo 
simulation.  For  example,  the  effect  of  varying  the  potential  barrier 

height  ( A1  mole  fraction)  is  shown  in  Fig.  2.20.  Although  the  peak 

velocity  increases  'with  the  barrier  height,  the  magnitude  of  the  negative 
differential  resistance  shows  a  maximum  for  a  barrier  height  of  about  0.2 
eV.  This  result  should  be  compared  with  the  one  obtained  from  the  electron 
temperature  model  (Fig.  2.7).  Fig.  2.21  shows  the  effect  of  changing  the 
layer  dimensions.  3y  increasing  the  ratio  of  the  AlGaA3  and  the  GaAs  layer 

width  the  peak  velocity  decreases  whereas  the  peak-to- valley  ratio 

increased.  The  larger  peak-to- valley  ratio  results  from  the  larger 
fraction  of  electrons  in  the  AlGaA3  as  shown  in  Fig.  2.22.  Again, 
Fig.  2.21  should  be  compared  with  Fig.  2.10  from  the  electron  temperature 
model. 


So  far  we  have  not  mentioned  the  effects  of  L  valley  in  the  GaA3.  The 
L  valley  i3  known  to  be  located  approximately  0.33  eV  above  the  T 
minimum[34, 36] .  If  the  barrier  height,  AE,  is  close  to  this  enrgy,  it  is 
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Velocity-field  characteristics  for  real  space  transfer  device 
with  AlGaAs  thickness  as  a  parameter  [64  ]  . 
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Fig.  2.22.  Relative  occupancy  of  electrons  in  AlGaAs,  with  AlGaAs  thickness 
as  a  parameter. 


possible  for  some  of  the  hot  electrons  to  transfer  to  the  L  valley  of  0aA3 
rather  than  to  the  A10aA3  layer.  Thi3  is  illustrated  in  Fig.  2.23  which 
shows  the  relative  populations  of  the  T  and  L  valleys  in  the  0aA3  and 

AlOaAs  for  the  parameters  used  in  Fig.  2.16.  With  this  set  of  parameters 
(A3  =  0.2  eV)  approximately  20  *  of  the  electrons  are  transferred  to  the  L 
valley  at  3  kV/cm.  If  the  band  edge  of  the  AlOaAs  is  aligned  with  the 

location  of  the  L  valley  (AS  =  0.33  eV),  the  effect  is  expected  to  be 
larger.  The  result  of  this  calculation  is  shown  '  in  Fig.  2.24.  Now  more 
than  40  %  of  electrons  transfer  to  the  L  valley.  This  figure  also  shows 
that  the  transfer  to  the  L  valley  is  slightly  more  efficient  than  real 
space  transfer  to  the  A1 Oa A3  layer.  For  the  real  space  transfer  mechanism 
to  be  dominant,  the  barrier  height  should  be  lower  than  ~  0.25  eV.  The 

effect  of  the  L  valley  in  the  AlOaAs  is  expected  to  be  very  3mall,  since 
the  electrons  in  the  low  mobility  AlOaAs  will  not  heat  up  as  much. 

The  Monte  Carlo  simulations  also  provide  the  distribution  function  by 
setting  up  an  energy  histogram[33 ] .  Fig  2.25  3hows  the  distribution 
function  sampled  at  the  GaAs  side  of  the  interface.  At  lower  fields  („i 
kV/cm)  the  curve  is  almost  a  straight  line  indicating  that  the  distribution 
is  approximately  Maxwellian.  As  the  field  increases,  however,  a 
significant  structure  appears  above  0.2  eV  which  corresponds  to  the 
conduction  band  edge  of  AlOaAs.  Since  the  mobility  in  the  AlOaAs  layer  is 
much  lower,  the  AlOaAs  electrons  do  not  heat  up  as  much  as  those  in  the 
OaAs  layer.  These  cooler  electrons  can  be  transfered  back  to  the  GaAs 
layer  causing  a  bump  in  the  distribution  function.  It  is  also  noted  that 
the  effective  electron  temperature  (inverse  slope  of  the  distribution 
function)  is  smaller  at  higher  energies  due  to  these  cooler  electrons  from 
the  AlOaAs.  This  situation  i3  quite  analogous  to  that  in  the  Gunn  effect, 
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Fig.  2.23.  Relative  occupancy  of  electrons  in  the  T  and  L  valleys  of 
GaAs  and  AlGaAs  versus  field  for  the  same  parameters  as  in 
Fig.  2.16. 
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Fig.  2.24.  Relative  occupancy  of  electrons  in  the  T  and  L  valleys  of  GaAs 
and  AlGaAs  versus  field.  All  the  parameters  are  same  as  in 
Fig.  2.20,  except  AE  =  0.33  eV. 
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Fig.  2.25.  Distribution  functions  for  electrons  in  the  GaAs  at  the 
GaAs-AlGaAs  interface. 


where  cooler  electrons  from  the  low-mobility  L  valleys  cause  a  bump  in  the 
distribution  function  in  the  f  valley  (compare  Pig.  2.25  with  Fig.  3  of 
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Ref.  33).  Also,  the  degree  of  the  deviation  from  the  straight  line 
indicates  the  degree  of  validity  of  the  electron  temperature  model,  which 
is  discussed  in  more  detail  in  the  next  section. 

2.3*4  Comparisons  with  the  electron  temperature  model 

As  seen  in  Fig.  2.25  the  distribution  function  can  be  considerably 
non-Max we Ilian  at  high  electric  fields.  In  this  case  the  electron 
temerature  model,  which  assumes  a  single  carrier  temperature  in  the  CaAs , 
is  not  adequate.  Moreover,  the  electron  temperature  model  is  based  on  more 
assumptions  (such  as  the  existence  of  the  quasi-Fermi  level)  than  the  Monte 
Carlo  simulations.  The  most  fundamental  difference  between  the  two  models 
i3  that  the  Monte  Carlo  method  is  a  one  particle  simulation  without 
carrier-carrier  interaction  being  taken  into  account,  whereas  the  electron 
temperature  model  assumes  strong  carrier-carrier  interaction  which  yields 
the  Maxwellian  distribution.  This  results  in  a  different  power  flow  rate 
as  is  explained  below. 

In  quite  general  situations  the  diffusion  of  energy  can  be 
accomplished  in  two  different  ways.  The  first  is  by  actual  carier 
transport,  '//hen  a  carrier  moves,  it  takes  its  energy  (kinetic  energy)  with 
it,  thereby  reducing  the  total  energy  left  behind.  The  other  is  by  the 
thermal  diffusion[42] .  This  i3  due  to  the  electron  temperature  gradient 
and  is  analogous  to  conventional  thermal  diffusion  of  gas  particles  in  a 
temperature  gradient.  These  two  effects  can  be  easily  demonstrated  for 
one-dimensional  case  by  multiplying  the  Boltzmann  equation  (3q.  (2.3))  by 
the  energy  and  integrating.  Cne  then  obtains(42,d3"! 
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where  <  is  the  thermal  conductivity  of  the  electrons  given  by 


k  =  2nkBD.  (2.17) 

Eq.  (2.16)  is  the  usual  power  balance  equation  with  two  additional  terms 
(inside  the  brackets)  due  to  the  inhomogeneity  of  the  system.  The  first 
term  represents  the  thermal  diffusion,  and  the  second  term  the  carrier 
transport  effect.  In  the  electron  temperature  model  both  of  these  effects 
are  included.  In  the  Monte  Carlo  simulation,  however,  only  carrier 
transport  is  included,  since  thermal  diffusion  i3  a  many-particle  effect. 
From  these  arguments  we  can  speculate  that  the  power  flow  out  of  the  CaAs 
is  larger  in  the  electron  temperature  model  than  in  the  Monte  Carlo  method. 
This  speculation  is  born  out  by  the  comparison  in  Fig.  2.26,  which  shows 
the  velocity-field  characteristics  from  the  two  models  using  the  identical 
set  of  parameters.  The  Monte  Carlo  simulation  3hows  a  large  negative 
differential  resistance,  "whereas  the  electron  temperature  model  result 
doesn't  even  show  velocity  saturation.  In  reality  the  characteristic  is 
expected  to  be  intermediate  between  the  two,  which  one  is  closer  depends  on 
the  carrier  density  in  the  system. 


2 .4  Two-dimensional  Effects 

2.1. 1  Size  quantization 


In  both  the  electron  temperature  model  and  the  Monte  Carlo  simulation 
we  have  considered  layer  dimensions  of  typically  400  .§  or  larger  to  avoid 
the  complexity  arising  from  two-dimensional  effects.  When  the  layer 
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dimension  i3  comparable  to  the  DeBroglie  wavelength  of  a  carrier,  bound 
3tates  are  formed  and  the  energy  levels  become  discrete  due  to  the 
restriction  of  the  carrier  motion  in  the  direction  normal  to  the  layer 
(referred  to  as  the  z-direction  in  this  section).  The  energy  levels  (for  a 
infinitely  deep  well)  are  then  given  by[44] 
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where  L  is  the  layer  thickness.  Each  value  of  S  is  connected  to  a 
continuum  of  levels,  called  a  subband.  In  GaAs-AlGaAs  heterostructure  this 
quantum  size  effect  (QSE)  is  expected  to  appear  for  layer  dimensions  below 
500  $[44].  The  quantization  of  carrier  motion  can  also  occur  at  the 
GaAs-AlGaAs  interface  due  to  band  bending[45 , 46 ] .  Electrons  or  holes  in 
such  a  system  can  be  considered  as  two-dimensional.  The  modification  of 
scattering  rate  due  to  this  two-dimensionality  and  its  influence  on  the 
real  space  transfer  mechanism  are  discussed  in  the  next  two  sections. 

2.4.2  Polar  optical  scattering  in  two  dimensions 

Polar  optical  scattering  is  the  most  dominant  scattering  mechanism  at 
low  energies  in  GaAs.  We  show  here  how  the  two-dimensional  effects  modify 
the  scattering  rate  of  thi3  mechanism.  Other  mechanisms  have  been 
considered  by  Hess[23]«  There  is  also  a  considerable  amount  of  work  on 


two-dimensional  scattering  in  connection  -with  silicon  MOS  inversion 
layers[47-49] • 


For  polar  optical  scattering  we  assume  that  the  phonon  modes  are 


three-dimensional,  since  the  OaAs  and  Al^a^  As  are  chemically  similar, 
especially  for  a  small  value  of  x.  We  then  start  with  the  Hamil tonian^oO] 


H  =  7  -  (  a  eiq‘r-  aT  e"lq,rl 
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where  a  i3  a  three-dimensional  phonon  vector,  a-*  and  a"*"  are  boson 

q  q 


operators,  and  C  is  given  by 


rfc2  , 

•  “°q  1 1  • 

1  2V  U  “  s  i 

—  oo  n  -1 


(2.21) 


We  only  consider  the  ground  state  for  the  electron  whose  vavefunction  is 


given  by 


(2.22) 


where  k„  and  r„  denote  the  components  parallel  to  the  interface,  L  is  the 
layer  thickness,  and  A  is  the  interface  area.  We  can  then  calculate  the 
matrix  element  of  the  electron- phonon  interaction: 
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Assuming  q  L  <  1  ,  we  obtain,  for  the  phonon  emission, 


(2.23) 
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The  total  scattering  rate  is  calculated  as 
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Performing  the  q^  integration  with  q^  <<  k„ ,  we  get[23] 
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where  upper  and  lower  eigne  correepond  to  phonon  absorption  and  emission, 
respectively,  and  the  notation  has  been  changed  (  k„  k,  and  q/;  q  )  .  To 
find  an  explicit  expression  for  x,  we  perform  the  integration  over  the 
3 -functional  ],  and  obtain: 
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for  phonon  absorption,  and 
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for  .phonon  emission.  In  these  expressions  3q  is  the  effective  field 
strength  for'  GaAs  (see  Append.  1),  and  X(x)  is  the  complete  elliptic 
integral  of  the  first  kind. 


In  Pig.  2.27  these  results  are  compared  with  the  polar  optical 
scattering  rate  in  three  dimensions.  As  seen  from  the  figure,  the 
two-dimensional  scattering  rates  are  considerably  higher  (for  both  phonon 
emission  and  absorption)  than  in  the  three-dimensional  case  for  any 
electron  energies.  An  enhancement  of  polar  optical  scattering  rate'  has 
been  experimentally  verified  by  Holonyak  et  al.[52]  using  photopumped 
multiple  quantum-well  GaA3-AlGaA3  heterostructures. 

2.4-3  Influence  on  the  real  space  transfer  mechanism 

The  enhancement  of  the  polar  optical  scattering  rate  alone  is  expected 
to  increase  the  threshold  electric  field  for  NDR  onset  by  a  factor  of  2  at 
room  tempera ture[5 ] •  It  is  rather  difficult,  however,  to  assess  all  the 
effects  of  two-dimensionality  on  the  real  space  transfer  mechanism  because 
of  it3  complexity.  For  a  realistic  device,  we  should  take  into  account  the 
existence  of  3ubbands.  The  energy  levels  of  these  subbands  can  be 
calculated  using  a  variational  principle[53> 54]  or  a  self-consistent 
methodr_53].  The  intervalley  and  the  polar  optical  scattering  rates  must  be 


obtained  including  those  subbands.  It  is  then  possible  to  calculate  the 
transport  property  in  the  two-dimensional  system  using  the  electron 
temperature  model[55],  or  the  Monte  Carlo  technique[o6].  The  transition 
from  two  to  three  dimensions  is  necessary  at  the  electron  energy  equal  to 
the  band  edge  of  AlCa As  to  account  for  the  real  space  transfer  mechanism. 
The  quantization  at  the  GaAs-AlGaAs  interface  due  to  band  bending[46]  is 
irrelevant  to  the  real  space  transfer,  because  when  the  carriers  become  hot 
they  are  no  longer  bound  at  the  interface. 

2 .5  Miscellaneous  Effects 

2.5.1  Band  bending 

3and  bending  effects  have  been  neglected  both  in  the  electron 
temperature  model  and  the  Monte  Carlo  simulation  for  the  sake  of 
simplicity,  although  it  is  possible  to  include  these  effects  in  the  Monte 
Carlo  simulation  in  the  form  of  a  position-dependent  electric  field.  The 
actual  potential  profile  including  band  bending  is  schematically  shown  in 
Fig.  2.23.  The  electric  field  created  by  the  ionized  donors  tends  to 
attract  energetic  electrons  and  pull  them  from  the  GaAs  into  the  AlGaAs. 
Hence,  it  is  expected  to  enhance  transfer  out  of  the  well  and  impede 

l 

transfer  back  into  the  well.  An  exact  treatment  to  calculate  the  band 
bending  should  employ  a  self-consistent  method[53]  to  account  for  the 
rearrangement  of  the  electrons.  However,  a  rough  estimate  can  be  given 
rather  easily  if  we  neglect  two-dimensional  effects  and  assume  a  uniform 
carrier  distribution  (valid  only  for  small  band  bending). 


Hi] 


Using  the  depletion  approximation!)  9 ] ,  we  obtain  ( Fig .  2.23) 
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where  n  is  the  carrier  density  in  the  GaAs.  Here  we  have  used  the  relation 
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(2.31 


at  the  start  of  the  transfer.  In  order  to  have  all  the  electrons  in  the 
3a As  layer  before  the  real  3pace  transfer,  we  want  to  have 
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Using  Sqs.  (2.29)  and  (2.90),  and  n  =  N^exp^  -.  ~uC )  ,  we  can  rewrite  thi3 


inequality  as 
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(2.33) 


If  we  replace  the  inequality  by  the  equality,  we  can  find  the  maximum 
allowable  layer  dimension  to  have  all  the  electrons  in  the  GaA3. 

Assuming  L^/L  1  =5,  =  1  x  10^  cm~^  and  AS  =  0.25  ev ,  we  obtain 

*  240  8,  h  “  0.12  M  ,  and  (Ji2  s  25  meV.  In  order  to  have  low-mobility 
AlGaAs,  the  layer  must  be  doped  higher  than  1  x  10l0  cm-^ .  Then  the  band 
bending  effect  is  certainly  not  negligible.  By  utilizing  compensation 


doping  of  the  Alla  As,  the  mobility  can  be  lowered  -without  large  band 
bending  effects.  In  any  case,  to  obtain  a  more  realistic  picture  of  the 
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rsal  space  transfer  mechanism,  band  bending  effects  should  be  examined  in 
the  Monte  Carlo  simulation.  Tunneling  through  the  barrier  may  also  be 
significant  for  large  band  bending.  Furthermore,  in  stead  of  trying  to 
eliminate  the  band  bending  effects,  there  might  be  a  possibility  of  charge 
storage  device  utilizing  this  effect. 


2.5.2  Transmission  coefficient  at  the  interface 


It  is  well  known[57,58]  that  when  an  electron  collides  with  a 
potential  step,  either  reflection  or  transmission  may  occur  with  finite 
probability.  Therefore,  at  the  GaAs-AlGaA3  interface  some  electrons  are 
transmitted,  ’while  others  are  reflected  during  the  transfer.  Wu  and 
Yang[59]  have  calculated  the  transmission  coefficient  across  heterojunction 
interfaces  taking  into  account  the  difference  of  effective  masses  on  the 
two  sides  of  the  junction.  Their  treatment  is  not  correct,  however, 
because  the  effective  mass  concept  is  not  valid  at  an  abrupt  potential 
step.  A  more  rigorous  calculation  has  been  done  by  Osbourn  and  Smith[oO] 
using  the  empirical  tight-binding  approximation.  Their  result,  as  compared 
with  Wu  and  Yang's  result,  is  3hown  in  Fig.  2.29.  The  figure  shows  the 
transmission  coefficient  at  the  interface  as  a  function  of  electron  energy 
for  a  wave  vector  perpendicular  to  the  interface.  Osbourn  and  Smith's 
result  in  terms  of  wave  vectors  has  been  replotted  using  the  relation 
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Their  calculation  gives  smaller  transmission  coefficient  than  obtained  by 
Wu  and  Yang's  calculation  as  pointed  out  by  Price[o1 ].  In  actual 
heteroj unctions  the  interface  i3  not  atomically  abrupt ,  but  is  graded  over 
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more  than  15  even  for  M3E  grown  heterojunctions[62].  This  compositional 
grading  has  been  shown  to  increase  the  transmission  coefficient 
3lightly[S5]. 

Glisson  et  al.[64]  have  compared  a  classical  transmission  model  and 
the  quantum-mechanical  transmission  using  the  formula  by  Wu  and  Yang  in  the 
Monte  Carlo  simulations.  In  the  classical  model,  the  transmission 
coefficient  is  a  step  function,  being  unity  at  the  conduction  band  edge  of 
the  AlGaAs.  Within  the  statistical  error  of  the  simulation  they  have 
observed  no  substantial  differences  in  the  results  between  these  two 
transmission  models.  If  the  results  by  Osbourn  and  Smith  are  used,  there 
may  be  a  small  difference.  Intuitively,  the  inclusion  of  the 
quantum-mechanical  transmission  is  nearly  equivalent  to  slightly  increasing 
the  barrier  height  in  the  classical  model.  In  any  case,  the  effect  is 
expected  to  be  small. 

2.5*3  Statistical  fluctuation 

Due  to  the  small  layer  dimensions,  the  statistical  fluctuations  of  the 

electron  density  and  the  impurity  density  might  be  substantial^^ ] •  For 

example,  if  we  assume  =  400  3,  =1  x  1013  cm-^ ,  and  the  device  area 

of  100  pm  by  100  pm,  the  total  number  of  electrons  will  be  N  =  4  x  10^ 

electrons,  which  is  rather  small  for  a  total  number  of  carriers.  The 

average  distance  between  electrons  in  this  case  i3  100  3.  For  the 

Impurity  density,  we  take  =  2000  8,  =  2  x  lO1"^  cm-"^ ,  and  the  same 
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device  area.  Then  the  total  number  of  impurities  in  the  AlGaAs  is  4  x  10. 
The  average  distance  between  impurities  is  170  S. 


Of  course,  this  problem  can  be  circumvented  by  utilizing  many  stacked 
layers.  By  the  "4BE  technique,  for  example,  it  is  possible  to  grow  more 


than  100  stacked  layers[44]. 


2 .6  Switching  Characteristics 


When  considering  the  application  of  the  NDR  mechanism  to  microwave  and 
switching  devices,  one  of  the  most  important  questions  to  ask  is  how  fast 
does  the  mechanism  occur  ,  or,  how  fast  does  it  respond  to  the  applied 
electric  field.  We  first  consider  the  transfer  time  of  the  electrons  from 
the  0aA3  into  the  AlGaAs  layer.  Under  collision- free  transport  conditions, 
the  transfer  is  due  to  the  thermionic  emission  current,  which  i3  given 
by[2S]  (see  Eq.  (2.12)) 
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where  n^  j.3  the  electron  'density  in  the  GaAs,  and  AE  is  the  potential 
barrier  height.  Comparing  Eq.  (2.55)  'with 
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gives  the  time  constant  t  as 
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Substituting  the  typical  values,  AE  =  0.2  eV,  kTg  =  0.17  eV,  and  L1  =  .100 
$,  we  get  ti_2=  4.8  x  10-1^  sec,  which  is  an  attractively  short  time. 
Since  the  time  required  for  heating  of  the  electrons  i3  much  larger  (~  5  x 


10-1^  sec)  than  ^ 
carriers. 


the  switch-on  time  is  limited  by  the  heating  of 


Depending  on  the  impurity  density  in  the  AlGaAs  layer,  the  transfer 
time  from  the  AlGaAs  to  the  GaAs  ^fter  switching  off  the  electric  field 
will  be  limited  either  by  diffusion  (high  density)  or  by  thermionic 
emission  (low  density).  For  the  diffusion  limited  case,  the  situation  is 
similar  to  charge  transfer  in  charge- coupled  devices,  where  electrons  move 
by  diffusion  from  one  gate  to  another  of  lower  potential  energy.  Using 
this  analogy,  we  obtain  the  transfer  time,T^  ,  from  the  AlGaAs  to  the 
GaAs,  given  by(5 , 66 ] 
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where  D  i3  the  diffusion  constant  in  the  AlGaAs.  This  formula  is  valid 
only  as  long  as  the  diffusion  concept  applies,  and  mean  free  path  for 
phonon  scattering ,  A^  >  i3  smaller  than  .  If  A  ^  is  longer  than  L1  ,  the 
probability  of  an  electron  being  captured  in  the  well  is  reduced  by  1^  A  ^ * 
'Then  we  have 
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?or  typical  values  such  as  1^  *  400  8,  *  2000  8,  A  ^  =  1000  8,  and  D 


u2kTL/e  =  10  Cm2/sec,  we  obtain  x2_1  =  4  x  1 0_1  1  sec.  This  give  an 


estimate  of  the  switch-off  time  for  the  diffusion  limited  case. 
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If  the  transport  is  collision- free,  the  time  T2-l  can  be  estimated  in  a 
similar  way  as  for  then  obtain 
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where  A  is  the  Richardson  constant  (120  A/cm  K  ),  m  is  the  effective  mass 
in  the  AlGaAs,  i3  the  free  electron  mass,  is  the  potential  created  by 
the  donors  in  the  AlGaAs  (  in  Fig.  2.25),  and  S  is  .the  effective  density 

4  C 

of  states  in  the  AlGaAs.  If  we  assume  eb/^T^  =  2,  we  obtain  T2_l  =  2.2  x 
1  P 

10  sec.  Therefore  the  switch-off  time  is  also  determined  by  tne  cooling 
time  of  electrons  in  the  GaAs. 


Of  course,  these  numbers  only  give  rough  estimates  of  the  transfer 
time.  A  more  rigorous  calculation  can  be  done  by  a  transient  Monte  Carlo 
method[o7]  including  the  real  space  transfer  mechanism.  It  is  also 
interesting  to  compare  the  real  space  transfer  time  and  the  transfer  to  the 
L  valleys,  and  observe  the  relative  population  as  a  function  of  time, 
especially  for  a  large  potential  barrier  height  (AS  -  0.3  eV). 

An  example  of  an  actual  device  configuration  for  a  switching  device  is 
shown  in  Fig.  2.30[68].  The  GaAs  and  the  AlGaAs  layers  are  contacted 
separately,  30  that  voltages  can  be  applied  independently  to  the  two 
materials.  Such  structures  can  be  realized  using  the  M3E  technique.  There 
will  be  a  negligible  current  in  the  AlGaAs  'in til  an  appreciable  fraction  of 
electrons  are  emitted  from  the  GaAs.  Fig.  2.31  illustrates  this  switching 
effect  calculated  by  the  electron  temperature  model.  The  electric  field  to 
the  AlGaA3  is  kept  constant  at  12  kV/cm.  The  GaAs  and  the  AlGaAs  currents 
can  be  switched  on  and  off  by  the  electric  field  applied  to  the  GaAs  layer. 


L iz01fjL  L2  =  1.0/x 
AE=250  meV 
/X2=  50cm2/V*sec 


Electric  Field,  F  (kV/cm) 


LP-1545 


Fig.  2.31.  Switching  characteristics  of  a  real-space  transfer  device.  The 
GaAs  and  the  AlGaAs  layers  are  contacted  separately.  The 
electric  field  to  the  AlGaAs  layers  is  kept  constant  at 
12  kV/cm.  Each  layer  can  be  switched  on  and  off  by  the  applied 
field  to  che  GaAs  layer. 
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The  dotted  curve  i3  obtained  in  the  case  of  no  thermal  conduction. 
Fig.  2.32  illustrates  another  aspect  of  the  switching  characteristics. 
Here  the  same  electric  field  is  applied  to  both  the  GaAs  and  the  AlGaAs 
layers.  The  figure  3hows  how  the  current  in  the  AlGaAs  increases  as  the 
field  increases  with  the  AlGaAs  mobility  as  a  parameter. 

The  possible  problems  which  might  be  encountered  in  the  application  of 
the  real  space  transfer  mechanism  to  switching  devices  include  the 
remaining  electrons  in  the  GaAs,  and  the  existence  of  traps  at  the 
interfaces  and  in  the  AlGaAs  layer.  Electrons  remaining  in  the  GaAs  layer 
above  threshold  field  cause  a  large  power  consumption.  Traps  cause  the 
loss  of  carriers  and/or  slower  switching  speed. 

2 .7  Summary 

In  this  chapter  we  have  described  the  mechanism  of  real  3pace  electron 
transfer  which  can  occur  in  GaAs-AlGaA3  heterostructures  .  A  negative 
differential  resistance  results  from  the  transfer  of  electrons  from  a 
high-mobility  GaAs  region  to  an  adjacent  low-mobility  AlGaAs  region  under 
the  influence  of  a  high  electric  field.  The  mechanism  is  analogous  in  many 
respects  to  the  Gunn  effect,  except  that  the  device  characteristics  can  be 
controlled  primarily  through  doping  densities,  layer  thicknesses,  and  the 
A1  mole  fraction  in  the  AlGaA3.  We  have  analyzed  the  mechanism  using  two 
different  methods,  the  electron  temperature  model  and  the  Monte  Carlo 
simulation.  Both  methods  clearly  illustrate  the  degree  of  control  of  the 
device  characteristics.  Some  comparisons  have  al30  been  made  between  the 
two  models.  The  electron  temperature  model  has  been  shown  to  give  a  larger 
power  flow  to  adjacent  layers  than  the  Monte  Carlo  simulation. 
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Two  dimensional  effects  might  be  important  in  actual  devices,  although 
they  have  been  neglected  in  the  two  methods.  More  work  needs  to  be  done  on 
this  subject.  Other  effects,  such  as  band  bending  and  quantum-mechanical 
transmission  have  been  discussed.  These  effects  might  also  influence  the 
actual  device  characteristics.  Especially,  band  bending  effects  should  be 
included  in  the  Monte  Carlo  simulations.  Switching  characteristics 
including  the  switching  speed  have  been  examined.  The  switching  speed  is 
estimated  to  be  on  the  order  of  1  x  10  sec,  which  is  an  attractively 

short  time  for  device  applications. 

There  is  also  a  possibility  of  enhancing  the  transfer-out  from  the 
0aA3,  such  as  simultaneous  application  of  an  electric  field  and  a  magnetic 
field  parallel  to  the  interface  to  divert  the  electron  trajectory.  Monte 
Carlo  simulations  including  a  magnetic  field  effect  should  be  easy.  The 
possibility  of  charge  storage  (resembling  charge-coupled  devices)  utilizing 
band  bending  effects  should  also  be  pursued.  Finally,  we  note  that  a 
similar  mechanism  is  important  in  other  semiconductor  devices-  For 
example,  in  silicon  MOS  devices  hot  electrons  move  to  the  Si-SiO^ 
interface,  causing  some  unwelcome  effects[70,71  ] .  In  the  same  way, 
electrons  are  emitted  into  substrate  in  FETs  causing  higher  output 
conduc tance[72 ] . 

Regarding  experimental  work  on  the  real  space  transfer  mechanism, 
Keever  et  al.[69]  have  recently  carried  out  measurements  of  the 
current-voltage  characteristics  of  3a A3- A 13a A3  heterojunctions.  They  have 
observed  current  saturation  and  negative  differential  resistance  in 
qualitative  agreement  with  the  Monte  Carlo  simulation  discussed  here. 
Their  results  tend  to  support  the  general  concept  of  real  space  transfer 
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described  in  thi3  chapter.  In  any  case,  there  is  no  question  that  hot 
electron  effects  in  semiconductor  heterostructures  will  attract  more 
interest  in  the  future. 


3.  BAND  STRUCTURE  DEPENDENT  TRANSPORT 
AND  IMPACT  IONIZATION  IN  la  As 

3 . 1  Introductory  Remarks 

A  large  number  of  semiconductor  devices  depend  on  high-energetic  (hot'' 
electrons  for  their  operation.  Impact  ionization  i3  an  essential  mechanism 
in  the  operation  of  many  semiconductor  devices,  such  as  photodetectors  and 
IMPATT  diodes.  At  present,  however,  the  theoretical  understanding  of  this 
effect  is  provided  by  a  number  of  theories[74-79]  which  contain  several 
adjustable  parameters  whose  physical  significance  is  not  well  known . 
Currently  the  most  widely  used  theory  of  impact  ionization  has  been  given 
by  3araff[76].  The  adjustable  parameters  of  hi3  theory  are  the  threshold 
energy  for  ionization,  the  optical  phonon  energy,  the  ionization  mean  free 
path,  and  the  mean  free  path  for  optical  phonon  scattering.  Although  some 
attempts  have  previously  been  made  to  determine  these  parameters 
theoretically[30,S1  ] ,  a  "complete'*  theory  of  Impact  ionization,  which  is 
capable  of  calculating  these  quantities  (and  therefore  the  ionization  rate) 
from  first  principles,  has  not  been  developed.  The  main  reason  is  that  any 
theory  applicable  at  high  electric  field  (causing  ionization)  must  abandon 
the  effective  mass  approximation  or  simple  extensions  using 
non-parabolicity  constant,  and  instead  include  a  realistic  band  structure. 
The  surprising  success  of  3araff's  theory  in  explaining  the  electric  field 
dependence  of  the  ionization  rate  is  due  to  the  adjustable  parameters, 
which  can  smear  out  the  band  structure  effects.  The  inclusion  of  the  band 
structure  in  analytical  solutions  of  the  Boltzmann  equation,  however,  is 
Impractical.  It  is  also  difficult  to  include  realistic  scattering 
mechanisms.  7or  example,  the  inclusion  of  both  small  angle  scattering  and 
randomizing  scattering  mechanisms  in  the  3ame  analytical  framework  is 
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difficult,  not  to  mention  several  different  mechanisms.  As  a  consequence 
all  the  previous  theories  are  only  applicable  to  specific  materials.  For 
instance,  Dumke's  theory[78]  is  only  applicable  to  InSb  or  InA3,  whereas 
Baraff's[76]  or  Chwang's[79]  treatment  is  only  valid  for  nonpolar 
materials,  such  as  Si  or  Ge. 

As  we  have  seen  in  the  analysis  of  real  space  transfer,  the  Monte 
Carlo  method  provides  an  alternative  to  the  solution  of  the  Boltzmann 
equation.  The  method  can  take  into  account  a  large  variety  of  scattering 
mechanisms,  and  therefore,  is  applicable  to  both  polar  and  nonpolar 
semiconductors.  It  can  calculate  the  quantities  of  interest  such  as  drift 
velocity,  mean  free  path,  average  electron  energy  without  any  a  priori 
assumptions  on  the  form  of  the  distribution  function.  Some  attempts  have 
previously  been  made  to  calculate  the  Impact  ionization  rate  by  the  Monte 
Carlo  method[33,84],  but  without  including  a  realistic  band  structure.  In 
this  chapter  we  describe  a  Monte  Carlo  method  that  includes  a  realistic 
band  structure  as  calculated  by  the  empirical  pseudopotential  method[35]* 
This  method  should  provide  a  tool  for  understanding  electronic  transport  in 
very  high  fields,  which  has  not  been  possible  with  conventional  methods 
utilizing  effective  mass  theory.  Of  course,  we  do  not  yet  have  all  the 
necessary  information  to  perform  a  rigorous  calculation  at  such  high 
fields.  For  example,  not  much  is  known  about  the  selection  rules  for 
scatterings  at  the  non- symmetric  points,  the  changes  of  the  ionization 
matrix,  and  the  scattering  rate  at  high  energies.  3ecause  of  this  lack  of 
information,  we  must  still  use  a  simpler  model  than  is  possible  with  this 
technique .  These  simplications  should  not  be  considered  as  restrictions  of 
the  method  itself.  As  more  information  becomes  available  in  the  future, 
the  method  can  be  improved  with  ease  to  accomodate  new  information.  In 


spite  of  the  simplications ,  the  results  we  obtain  provide  insight  into  how 
the  various  parameters  in  previous  theories  are  connected  and  into  how  the 
band  structure  influences  the  impact  ionization  rate.  In  particular,  the 
method  is  applied  to  study  the  electron  initiated  ionization  in  GaAs. 
Finally  the  results  also  give  us  information  about  the  accuracy  of  the 
pseudopotential  band  structure  at  high  energies. 

3 .2  Summary  of  Experimental  Results 

Before  we  discuss  the  theory  of  Impact  ionization,  it  is  instructive 
to  summarize  the  available  experimental  data.  Various  experimental 
techniques  to  measure  the  ionization  rate  are  described  in  detail  in  a 
review  paper  by  Stillman  and  Wolfe[3].  Their  article  also  contains  some  of 
the  experimental  data  on  the  electron  initiated  ionization  rate  in  GaAs. 
Fig.  3-1  summarizes  more  recent  data[36-90].  The  experimental  data  usually 
shows  a  1/E  or  1  /E2  dependence  of  the  ionization  rate.  As  can  be  seen  from 
the  figure,  the  data  of-  different  workers  scatter  almost  by  an  order  of 
magnitude . 

Of  special  interest  are  the  results  of  Pearsall  et  al._39],  who 
measured  the  electron  ionization  rate  with  the  electric  field  applied  in 
three  different  crystallographic  directions.  Their  data  are  replotted  in 
Fig.  3.2.  They  have  measured  the  highest  ionization  rate  in  the  <110> 
direction  and  the  lowest  in  the  <11  1>  direction.  They  have  attributed  this 
difference  to  ballistic  electrons  and  electron  tunneling  to  upper 
conduction  band[39,91  ].  Although  these  data  raise  an  interesting  question 
as  to  how  the  band  structure  actually  influences  the  ionisation  rate,  their 
notion  of  ballistic  electrons  3eems  to  be  incorrect,  as  is  shown  in  this 
work.  More  systematic  and  reliable  data  are  necessary  to  make  a  comparison 
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Fig.  3.1.  Experimental  ionization  rates  of  electrons  for  GaAs  at  room 
temperature  [86-90  ]  . 
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with  theory. 

3 .3  Previous  Theories  of  Impact  Ionization 

Wolff[74]  was  the  first  to  calculate  the  ionization  rate  in 
semiconductors.  ‘He  applied  the  gas  discharge  theory  to  solve  the  Boltzmann 
equation  taking  into  account  the  effect  of  electron- phonon  and 
pair- producing  collisions  on  the  distribution  functions.  The  velocity 
distribution  function  was  approximated  as 

n  (v,9)  *  nQ(v)  +  n^(v)  cos0  ,  (3.1) 

where  v  is  the  electron  velocity,  and  9  is  the  angle  between  the  velocity 
and  the  electric  field.  This  is  an  energy  diffusion  theory,  in  which  the 
electrons  undergo  many  collisions  when  moving  to  higher  energies.  The 
Boltzmann  equation  was  then  solved  to  calculate  the  ionization  rate  with 
the  result 

ct(  F)  ~  exp(-A/'F^)  ,  (3.2) 

where  F  is  the  electric  field. 

Shockley[75] ,  on  the  other  hand,  argued  that  ionization  is  mainly  due 
to  "lucky"  electrons  which  completely  escape  phonon  scatterings  and  reach 
the  threshold  energy.  In  this  streaming  approximation  the  distribution  is 
a  3pike  in  the  direction  of  the  electric  field.  He  considered  the  relative 
probability  of  phonon  scattering  and  pair  production,  and  obtained  an 
ionization  rate  whose  dependence  on  F  is  given  by: 

a(  F)  ~  exp(-B/'F)  .  (3.3) 
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Next,  Baraff[76 , 92  ]  solved  the  time  dependent  Boltzmann  equation  and 
showed  that  his  result  contained  Shockley's  result  as  a  low  field  limit,  and 
Wolff's  result  as  a  high  field  limit.  His  theory  gives  the  "universal" 
curves  with  phonon  mean  free  path  and  ionization  threshold  energy  as 
parameters,  which  are  adjusted  to  fit  theory  to  experimental  data. 
However,  it  does  not  provide  a  way  to  calculate  these  quantities,  nor  does 
it  include  the  band  structure. 

Recently,  Chwang  et  al.[79]  took  a  different  and  interesting  approach 
using  a  finite  Markov  chain  formation.  Their  method  is  based  on  the 
calculation  of  a  transition  matrix  which  characterizes  the  transition 
probability  between  virtual  states  defined  by  small  discrete  ever gy 
intervals.  Interesting  as  it  is,  their  method  is  unfortunately  limited  by 
the  analytical  formulations.  It  still  requires  the  same  assumptions  as 
Baraff's  theory  and  does  not  produce  much  more  information.  For  example, 
an  assumption  of  a  constant  mean  free  path  for  phonon  scattering  is  3till 
necessary.  Moreover,  the  Markov  formulation  is  only  applicable  to  nonpolar 
semiconduc  tors . 

Nevertheless,  Baraff's  and  Chwang' s  theories  contain  some  "truth" 
about  the  impact  ionization  mechanism,  as  does  Shockley's  or  Wolff's 
approach.  How  they  are  related,  and  how  they  complement  each  other  will  be 
clear  as  a  result  of  the  Monte  Carlo  calculation  described  in  this  work. 
This  Monte  Carlo  method  includes  a  realistic  band  structure.  As  a  result, 
the  orientation  dependence  of  the  Impact  ionization  rate  can  be  calculated 
for  the  first  time. 
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3  «4  The  Band  Structure  of  GaAs 

The  band  structure  of  GaAs  has  been  calculated  using  the  empirical 
pseudopotential  method  as  described  by  Cohen  and  Bergstresser[35] •  Only 
the  lowest  conduction  band  has  been  considered.  The  effect  of  higher  band 
is  briefly  discussed  in  the  later  sections.  Advantage  is  taken  of  the 
43- fold  symmetry  of  the  Brillouin  zone  of  the  zinc-blend  structure[93 ].  It 
is  then  necessary  to  examine  only  that  1/48  th  part  of  the  zone  indicated 
in  Fig.  3* 3>  This  region  is  defined  by  the  conditions: 

0  <  k  <  k  <  k  <1,  (3.4) 

—  z  —  x  —  y  — 

and 

k  +  k  +  k  <  ,  (3.5) 

x  y  z  —  2 

where  all  the  k  components  are  in  units  of  2Va  (a  is  the  lattice  constant; 

a  -  5-64  8  for  GaAs).  Mesh  points  (k  k  ,  k  =  0.0,  0.1,  ..)  are  sampled 

x  y  z 

from  this  region,  and  the  energy  and  its  gradient  (velocity)  at  each 
k-point  is  calculated.  A  total  of  249  points  have  been  sampled,  with  156 
points  within  wie  sampling  region.  The  extra  93  points  outside  the  region 
are  necessary  for  the  interpolation  of  energy  in  the  proximity  of  the 
surface  of  the  sampling  region.  Table  3.1  lists  the  calculated  S(k) 
relation  for  these  249  points.  Fig.  3-5  illustrates  the  isoenergy  Tines  in 
the  cross  section  of  the  Brillouin  zone  shown  in  Fig.  3*4,  with  the  numbers 
representing  the  electron  energy  from  the  bottom  of  the  conduction  band  (T 
point).  It  can  be  seen  that  the  f  valley  is  nearly  isotropic,  whereas  the 
X  valleys  are  more  elliptic.  A  similar  plot  of  isoenergy  lines  in  the 
cross  section  shown  in  Fig.  3-6  is  illustrated  in  Fig.  3*7.  The  overall 
electron  energy  in  this  cross  section  is  considerably  lower  than  in  the 


Table  3.1  (a).  E(k)  relation  for  the  lowest  conduction  band  of  GaAs 
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5.  Isoenergy  lines  of  the  lowest  conduction  band  of  GaAs  in  the 
cross  section  shown  in  Fig.  3.--.  The  numbers  represent  the 
energies  measured  from  the  T  minimum. 
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Fig.  3.7.  Isoenergy  lines  of  the  lowest  conduction  band  of  GaAs  in  the 
cross  section  shown  in  Fig.  3.6.  The  numbers  represent  the 
energies  measured  from  the  T  minimum. 
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cross  section  of  Fig.  3-4.  This  band  structure  i3  stored  in  memory,  and 
used  in  the  Monte  Carlo  simulation  which  is  described  in  Section  3.6. 


In  the  study  of  ballistic  electron  transport  in  three  major  crystal 
directions,  i.e.,  the  <100>,  <110>,  and  < 1 1 1 >  directions,  extra  k-points 
have  been  sampled,  and  the  S(k)  relation  is  calculated.  The  result  of  the 
calculation  is  shown  in  Fig.  3*8.  It  is  obvious  from  the  figure  that  the 
use  of  effective  mass  and  non-parabolicity  constant  is  not  valid  for 
electron  energies  above  approximately  1  eV  in  some  directions.  In  fact, 
the  effective  mass  defined  as 


(3.6) 


_1_  =  J_  3  E(k) 
■k  1 2  ,,  2 

m  n  3k 


goes  to  negative  values  at  higher  energies. 


3-5  Ballistic  Electron  Transport  and  Phonon  Seal 


The  term  "ballistic  electrons"  is  used  to  denote  those  electrons  which 
do  not  suffer  phonon  scattering.  This  is  equivalent  to  the  "lucky" 
electron  notion  in  Shockley's  theory.  Since  the  possible  contribution  of 
ballistic  electrons  to  impact  ionization  has  been  suggested[91  ] ,  the 
behavior  of  ballistic  electrons  has  been  examined  using  the  pseudopotential 
band  3tructrure[94].  The  study  has  been  performed  by  solving  the  equations 


h  =  e? 
dt 


(3.7) 


90 


v  =  £  7kE(^)  (3.8) 

-►  -*■ 

where  ?  is  the  applied  electric  field,  k  is  the  electron  wave  vector,  E  is 
the  electron  energy,  and  v  is  the  group  velocity  of  the  electron.  With  the 
initial  condition  k  =  0  at  t  ■  0,  we  obtain 

hk  =  et?  .  (3.9) 

The  field  is  assumed  to  be  constant.  Eqs.  (3*3)  and  (3-9)  are  solved 
simultaneously  to  express  v  and  E  as  a  function  of  time  t.  The  results  of 
the  calculations  are  shown  in  Pigs.  3*9  and  3*10  for  the  three  major 
crystallographic  directions.  The  electric  field  ha3  been  chosen  to  be  509 
kV/cm,  a  typical  field  for  impact  ionization.  Fig.  3*9  shows  the  electron 
velocity,  v,  as  a  function  of  time.  The  orientation  dependence  of  the 
ballistic  behavior  is  obvious  from  this  figure.  The  highest  peak  velocity 
is  reached  in  the  <100>  direction  (~  1.1  x  10®  cm/sec)  and  the  lowest  in 
the  < 1 1 1 >  direction  („  0.5  x  108  cm/sec).  Fig.  3-10  shows  the  variation  of 
electron  energy  with  time  as  measured  from  the  conduction  band  edge.  The 
rate  of  increase  i3  largest  in  the  <100>  direction  and  smallest  in  the 
<1 1 1 >  direction. 

i 

In  reality,  however,  ballistic  transport  must  compete  with  scattering 
processes.  It  will  be  shown  by  the  Monte  Carlo  simulation  that  on  the 
average  an  electron  can  travel  ballistically  for  only  ~  3  x  14  sec 
before  it  suffers  a  phonon  scattering.  In  the  <11 1>  direction,  the 
electron  can  never  gain  sufficient  energy  ballistically  for  impact 
ionization[39] .  In  the  <100>  direction,  the  impact  ionization  threshold 
can  be  reached  only  if  electrons  tunnel  in  k-space  to  the  next  higher  band 
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-  0.2  eV  above  the  principal  conduction  band[39].  Therefore,  scattering 
events  become  crucial  for  the  occurrence  of  impact  ionization  in  these 
directions . 

Electrons  can  be  scattered  to  other  regions  of  the  3rillouin  zone  with 

a  single  scattering  event  being  sufficient  to  permit  the  electrons  to  reach 

threshold  energy.  This  mechanism  i3  illustrated  in  Figs.  3.11  and  3-12  for 

an  electric  field  applied  in  the  < 1 1 1 >  direction.  An  electron  starts  at 

the  r  point  and  moves  along  the  <  1  1 1  >  direction.  At  point  A  (k  - 

•^(0. 3, 0.3,0. 3))  i  the  energy  is  at  the  maximum  for  this  direction,  but  it  i3 
a 

still  much  less  than  the  threshold  energy.  Subsequently ,  the  electron  can 
be  scattered  (by  a  phonon  or  impurity)  to  some  other  point  in  the  Brillouin 

zone,  point  B,  for  example.  Following  this  scattering  event,  the  <11 1> 

component  of  the  electron  wave  vector  continues  to  increase.  However,  the 
wave  vector  points  in  a  direction  different  from  the  <11 1>  so  that  the 

electron  can  now  reach  a  higher  energy.  As  shown  in  Fig.  3*12,  the 

electron  can  actually  exceed  the  threshold  energy  for  Impact  ionization  (-. 
2.0  eV). 

This  i3,  of  course,  only  one  example  of _an  electron  trajectory  to  show 

the  Importance  of  scattering  processes  to  impact  ionization.  The  actual 

calculation  of  the  impact  ionization  rate  must  involve  averaging  of  all  the 

possible  electron  trajectories  until  the  electron  reaches  the  threshold 

* 

energy.  Thi3  is  achieved  by  the  Monte  Carlo  method  which  is  described  in 


the  next  section. 


(a)  (110)  section  of  the  first  Brillouin  zone. 

(b)  Wave  vector  trajectory  of  electron  in  this  plane  under  the 
influence  of  electric  field  in  the  <111>  direction. 
Electron  is  scattered  from  A  to  B.  Energy  change  in  the 
scattering  process  has  been  neglected. 


3*6  Simulation  Method 
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The  Monte  Carlo  simulation  keeps  track  of  an  electron  k-vector  in  the 
Brillouin  zone  until  it  reaches  the  threshold  energy  for  impact  ionization. 
This  is  done  with  a  knowledge  of  scattering  mechanisms,  scattering  rates, 
and  band  structure  in  the  whole  Brillouin  zone.  The  3(k)  relation  for  an 
arbitrary  k-point  can  be  calculated  in  the  following  way.  First,  the 
k-point  is  mapped  into  the  sampling  region  (Fig.  3-3)  by  symmetry 
operations  which  consist  of  taking  absolute  signs  and  the  permutation  of 
wave  vector  components.  The  energy  is  then  calculated  by  quadratic 
interpolation  utilizing  the  energies  and  the  gradients  of  the  surrounding  3 
mesh  points.  The  gradient  is  interpolated  only  linearly.  Applying  the 
inverse  operations  on  the  calculated  energy  and  gradient  gives  the  3(k) 
relation  and  the  gradient  at  the  original  k-point. 

Next  we  need  to  know  the  phonon  scattering  rate.  Ideally  the 
scattering  rate  should  be  calculated  at  each  k-point  in  o-der  to  take  into 
account  the  overlap  integral[95]  (see  Sq.  ( A 1 . 7 )  in  Append,  l).  Also,  when 
the  initial  or  final  electron  state  is  not  on  the  symmetry  points,  the 
selection  rulesl_96]  become  less  restrictive  and  may  give  rise  to  additional 
scattering  mechanisms.  Moreover,  even  near  the  bottom  of  the  valleys,  it 
is  known  that  the  scattering  rates  are  different  in  the  "  ,  L,  and  X 
valleys.  In  spite  of  these  facts,  we  have  assumed  the  scattering  rate  to 
be  isotropic  (only  energy-dependent)  for  simplicity  and  because  of  lack  of 
additional  information.  'Ve  have  taken  the  scattering  rate  as  given  for  the 
central  valley.  Thi3  overestimates  the  scattering  rate  when  the  electron 
is  in  the  satellite  valleys.  The  simplification  is  partly  .justified  by  the 
fact  that  the  scattering  rate  of  different  valleys  approach  each  other  at 


I 
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higher  energies.  Furthermore,  3inoe  scattering  to  upper  bands  is  possible 
in  reality  (which  increases  the  scattering  rate  in  the  satellite  ’/alleys'1, 
the  overestinacion  is  at  least  partly  compensated. 

The  values  of  the  parameters  for  the  calculation  of  the  scattering 
rate  are  the  same  as  the  ones  used  in  the  simulation  of  the  real  space 
transfer  mechanism,  which  are  known  to  give  a  good  fit  to  experimental  da:a 
of  the  Gunn  effect[34-].  Below  0.33  eV  only  polar  optical  scattering  occurs 
in  the  central  valley.  Above  0.33  eV  polar  optical  scattering  occurs  only 
when  an  electron  is  in  the  central  valley  defined  arbitrarily  as 

-  0.3  <  k  ,  k  ,  k  <  0.3  ,  (3.10) 

—  x  y  z  — 


where  the  components  are  in  units  of  Tm/a.  Otherwise  intervalley 

scattering  occurs.  It  is  not  appropriate  to  simply  extend  the  scattering 

rate  to  higher  energies  because  of  the  complicated  band  structure.  Because 

the  intervalley  scattering  rate  is  proportional  to  the  density  of  final 

states,  and  the  density  of  states  in  the  conduction  band  decreases  nearly 

quad ratically  above  1.5  eV[97],  we  have  assumed  a  quad ratically  decreasing 

scattering  rate  above  1.5  eV.  The  resultant  total  scattering  rate  as  a 

function  of  electron  energy  is  shown  in  Fig.  3.13  (solid  line''.  The 

1 4  -1 

maximum  scattering  rate  is  4*5  x  10  sec  at  1.5  eV. 


In  the  <100>  direction  the  threshold  state  for  electron  initiated 
ionization  lies  in  the  second  conduction  oand^BS].  An  electron  can  funnel 
through  the  "pseudo-gap"  (~  0.2  eV)  between  the  lowest  and  the  second 


conduction  band  to  reach  threshold.  Mo  attempt  has 
this  tunneling  mechanism.  Since  the  tunnelinsr  time 
the  order  of  1  x  1 0_1 ^  sac[39,93],  and  the  inter/ all 
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an  electron  energy  of  2-0  e7  is  much  shorter  (<  1  x  1  0_1  ^  sec)  than  this 
tunneling  time,  electrons  are  more  likely  to  be  3-c  arse  red  before  they  can 
tunnel  to  the  upper  band.  Therefore,  the  contribution  of  these  tunneling 
electrons  to  impact  ionization  is  expected  to  be  small. 

The  final  state  of  the  scattering  process  is  determined  in  the 
following  way.  Since  polar  optical  scattering  is  dominant  only  at  low 
energies  in  the  central  val'ley,  the  usual  formula  (see  Append.  0  with 
effective  mass  and  non-parabolicity  terms  is  used  to  choose  a  candidate  for 
the  final  k-point.  The  energy  at  this  k-point  is  then  recalculated  using 
the  exact  band  structure  to  check  if  it  is  within  an  allowed  range 
(typically  30  meV)  around  the  final  energy  (for  example,  3+"h“0  in  the  case 
of  phonon  absorption).  If  it  is  outside  this  range,  a  different  final 
state  is  chosen  and  the  process  is  repeated  ’until  a  proper  3tate  within  the 
correct  energy  range  is  found.  Intervalley  scattering  is  known  to  be 
completely  randomizing^ 2 6 ] .  For  this  mechanism,  once  the  final  energy  is 
calculated,  those  mesh  points  whose  energies  are  within  the  allowed  range 
are  tabulated.  One  of  them  is  then  randomly  selected  as  the  final  state. 

Our  treatment  of  phonon  scattering  processes  represents  a  compromise 
between  accuracy  and  numerical  tractability.  For  a  finite  number  of  mesh 
points,  the  energy  separation  between  any  two  k-points  is  finite.  For 
example,  for  our  156  mesh  points  this  energy  separation  can  be  as  large  as 
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final  energy  can  be  infinitely  small.  The  number  of  k-points  in  this 
energy  range  for  a  given  final  energy  is  proportional  to  the  density  of 


states  at  each  region  of  k  space  with  thi3  final  energy.  Therefore,  the 
procedure  is  basically  correct,  since  the  scattering  rate  for  deformation 
potential  scattering  is  proportional  to  the  final  density  of  states. 


The  impact  ionizing  collision  is  treated  as  an  additional  scattering 
mechanism.  We  assume  an  isotropic  threshold  energy  of  2.0  e7.  Anderson 
and  Crowell[so]  have  shown  that  the  threshold  energy  actually  depends  on 
the  k-vector.  However,  their  graphical  procedure  is  almost  impossible  to 
perform  in  three-dimensional  momentum  space.  A  more  systematic  approach 
nay  be  possible]?? ] .  If  the  threshold  energy  is  calculated  for  each 
k-point,  it  can  be  easily  included  in  this  simulation  procedure.  The 
impact  ionization  probability  can  be  calculated  from  the  matrix  “"dement  for 
the  screened  Coulomb  interaction[34, 100, 101  ].  However,  here  we  use  a 
simpler  model  demonstrated  by  Xeldysh]77 , 1 02 ]  and  used  by  others1]??]. 
According  to  Keldysh  the  probability  of  impact  ionization  can  be 


represented  as 


t.(E) 


_j_  P 

T(Ei>  !-  Ei- 


(3.11) 


wnere  E  is  the  electron  energy,  Hf  is  the  threshold  energy,  ■/  x  '3^'  is  the 
scattering  rate  at  E  =  E^,  and  ?  i3  a  dimensionless  constant  which  i3 
usually  much  larger  than  'unity.  This  formula  is  valid  for  semiconductors 
with  large  dielectric  constants.  We  take  ?  as  a  parameter.  ?  =  50  -  » 
has  been  used  by  Chwang  et  al.]??].  As  shown  by  3araff]?S ' ,  and  then  by 
Chwang]??]  the  impact  ionization  rate  does  not  strongly  depend  on  this 
parameter  as  long  as  ?  is  large  compared  with  unity.  The  energy  dependence 


of  impact  ionization  probability  f o r  ?  = 


illustrated  in 


400  is 


101 


Once  the  scattering  rata  ani  the  ionization  probability  are 
determined,  the  re3t  of  the  simulation  procedure  is  similar  to  the  one 
described  for  the  real  space  transfer  mechanism,  except  that  instead  of 
determining  the  drift  time  by  the  integral  (3q.  (2.14)),  the  scattering 
probability  (At/x(E))  is  calculated  at  each  time  interval,  At,  and  compared 
with  a  random  number.  This  is  necessary  because  of  the  complicated  E(  k) 
relation.  At  is  taken  to  be  approximately  1/10  th  the  average  drift  time. 
The  simulation  starts  by  releasing  an  electron  with  zero  energy  at  the 
bottom  of  the  central  valley.  The  energy  and  the  k  vector  of  the  electron 
are  traced.  When  Impact  ionization  occurs,  the  energy  is  reinitialized  to 
zero  to  start  a  new  history.  This  i3  justified  by  the  fact  that  the 
resultant  electron  after  ionization  lies  very  close  to  the  bottom  of  the 
central  valley[30].  The  impact  ionization  rate  can  be  obtained  by 
averaging  each  distance  that  an  electron  travels  until  impact  ionization 
occurs  over  a  sufficient  number  of  ionizations.  The  distance,  Ax,  traveled 
during  each  drift  is  calculated  either  by  accumulating  a  differencial 
distance,  vAt,  or  by  utilizing  the  relation 

AE  =  eFAx  ,  (3.12) 

where  AS  is  the  energy  gained  during  the  drift.  The  velocity  v  is 
calculated  from  the  gradient  of  S(k)  relation. 

3.7  Results 


3.7.1  Contribution  of  ballistic  electrons 


3y  terminating  the  simulation 
suffers,  the  behavior  of  ballistic 
we  can  determine  the  extent  which 
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electrons  contribute  to 
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impact  ionization.  Since  there  is  no  electron  initiate!  tnresnold  3tate  in 
the  <100>  or  the  <11  1>  direction!.  89] ,  we  only  consider  the  <110>  direction. 
We  have  also  changed  the  threshold  energy  to  1  .7  eV,  which  is  the  correct 
threshold  energy  in  this  direction[ 8 0 ] .  Typically  100,000  trials  have  been 
done  for  each  electric  field. 

The  result  of  the  calculation  shows  that  an  electron  travels  on  the 
average  approximately  200  8,  for  an  average  time  of  3  w  1 0-'' ^  sec  before 
the  first  scattering  event.  These  numbers  differ  slightly  for  different 
orientations.  3 y  counting  those  electrons  which  cause  impact  ionization 
instead  of  scattering,  we  can  estimate  the  contribution  of  ballistic 
electrons  to  the  impact  ionization  rate.  If  the  same  scattering  rate  i3 
used  as  shown  in  Tig.  3-13  (known  to  give  a  good  fit  to  the  Ounn 
effect[34]),  we  find  no  electrons  (less  than  0.001  %)  causing  impact 
ionization.  There  may  be  some  uncertainties  in  the  scattering  rate, 
particularly  in  the  values  for  the  effective  masses  or  the  deformation 
potential  constants.  To  find  the  maximum  possible  contibution  of  ballistic 
electrons  to  Impact  ionization,  a  smaller  scattering  rate  has  been  tried. 
We  have  U3ed  the  values  given  by  Vinson  et  al/l03j5  =  0.4  eV; 

3  =  Q.38  eV ;  D r  T  =  1  . 1  x  1  09  eV/om;  and  D„  T  •=  2.S  x  1  03  eV/cm. 

Thi3  gives  the  scattering  rate  3hown  by  the  broken  line  in  Fig.  3*13.  This 
rate  is  approximately  half  of  the  previous  value.  Using  this  scattering 

rate  in  our  calculation  we  obtain  the  results  shown  in  Fig.  3.14.  This 

figure  3how3  the  probability  that  an  electron  causes  impact  ionizatior 
prior  to  its  scattering  by  a  phonon  as  a  function  of  electric  field.  As 
can  be  3een,  even  at  the  maximum  (300  kV/cm  for  P  =  00 )  ,  only  0.2  -1  of  th- 
electrons  causing  Impact  ionization  are  "ballistic".  Therefore  we  conelud- 
that  the  contribution  of  ballistic  electrons  to  impact  ionization  i: 


Probability  of  Ballistic  Ionization 
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negligibly  small  if  our  present  ’understanding  of  the  scattering  rate  is 
correct.  This  conclusion  negates  the  considerations  by  Dapasso  et  al.[9l] 
who  suggested  that  ballistic  electrons  give  a  non-negligible  contribution 
to  the  total  ionization  rate. 

These  "ballistic"  electrons  were  discussed  by  Shockley  who  called  them 
"lucky"  electrons[75 ] .  We  have  shown  that  Shockley's  theory  gives 
ionization  rate  that  are  too  small.  It  is  interesting  to  note,  however, 
that  the  two  curves  in  Fig.  3*14  show  the  correct  1 /E  dependence  as  in 
Shockley's  theory  inspite  of  the  more  complicated  band  structure  and 
scattering  rate  that  we  used. 

3-7.2  Transport  properties  and  ionization  rate 

In  the  calculation  of  the  impact  ionization  rate  a  typical  simulation 
consists  of  approximately  200,000  to  400,000  scattering  events.  Depending 
on  electric  field  this  would  give  40  -  300  impact  ionization  events. 
Fig.  3*15  shows  a  typical  trajectory  of  the  k  vector  in  the  Brillouin  zone 
for  an  electric  field  of  500  kV/cro  in  the  <100>  direction.  The  solid  lines 
represent  the  drift  of  the  electron,  and  the  broken  lines  represent  the 
scatterings  from  one  end  to  the  next.  When  the  k  vector  lies  outside  of 
the  Brillouin  zone,  it  is  placed  back  inside  the  zone  to  the  equivalent 
point.  This  is  done  by  adding  the  appropriate  reciprocal  lattice  vector  to 
the  original  k  vector.  As  seen  from  the  figure  the  drift  time  is  very 
short  because  of  the  high  scattering  rate  at  higher  energies.  The  electron 
is  frequently  scattered  over  practically  the  entire  Brillouin  zone. 
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Pig.  3*16  shows  how  the  calculated  drift  velocity  converges  as  the 

number  of  scatterings  is  increased.  Following  some  initial  fluctuations, 

the  velocity  seems  to  "settle  down"  after  approximately  2000  scatterings. 
Therefore,  200,000  scatterings  should  give  a  good  estimate  of  the  drift 
velocity.  Similar  curves  are  shown  in  Fig.  3*17  for  the  ionization  rate. 
Because  of  the  limited  CPU  time  it  has  not  been  possible  to  take  averages 
over  mors  than  300  ionizations.  However,  the  convergence  is  fairly  good 
after  10  ionizations.  From  thi3  figure  the  statistical  fluctuation  is 
estimated  to  be  appo roximatly  20  't.  The  problem  of  statistical  fluctuation 
can  be  overcome  by  repeating  the  simulation  only  for  the  high  energy 
tail[l04].  This  has  not  been  attempted  in  this  work. 

Fig.  3*19  shows  the  average  electron  energy  as  a  function  of  the 

electric  field.  The  reason  for  the  steeper  increase  beyond  IOC  kV/cm  is 
not  known,  but  may  be  related  to  the  band  structure.  Fig.  3-19  shows  the 
electron  mean  free  path  as  a  function  of  electric  field.  In  the  electric 
field  range  where  Impact  ionization  occurs,  the  mean  free  path  changes  from 
50  §  to  30  3.  This  is  in  good  agreement  -with  the  experimental 

iata[9,105]  and  the  conventional  Monte  Carlo  calculation[31  1 .  The  reason 
that  our  calculation  agrees  in  this  respect  'with  the  conventional  method 
which  does  not  include  a  realistic  band  structure  is  that  the  mean  free 
path  is  main1y  determined  by  the  average  electron  energy,  which  i3  still 
3saall  enough  ( ~  0.3  eV)  for  effective  mas3  and  non-parabolicity  corrections 


to  be  sufficient, 


The  calculated  electric  field  dependence  of  the  electron  drift 
velocity  is  shown  in  Pig.  3*20.  The  broken  curve  represents  the 
experimental  data  by  Ruch  and  Xino^l 06]  at  low  electric  field  (<  1  -  kV  ' rz) , 
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.  3.20.  Calculated  electron  drift  velocity  in  GaAs  at  room  temperature 
compared  with  the  experimental  data _J100 , 107  ]  .  The  calculated 
values  are  from  the  slope  of  the  E(k)  curve  (!)  or  from 
Eq.  (3.12)  (o). 


and  by  Houston  and  3vansM07j  at  high  field  {20  -  1 00  '<’// on).  The 
agreement  is  good  over  the  entire  range  of  electric  fields  experimentally 
investigated .  The  result  using  Fq.  (3.1l)  gives  much  better  fit  than  the 
result  using  the  slopes  of  the  3(*c)  relation.  It  is  suspected  that  the 
accumulation  of  numerical  errors  in  the  3lope  calculation  is  responsible 
for  the  discrepancy.  The  slight  deviations  between  theory  and  experiment  at 
higher  fields  are  believed  to  be  mainly  due  to  the  pseudopotential  band 
structure  which  gives  the  satellite  valley  effective  massses  larger  than 
are  usually  measured.  In  Fig.  3-20  it  can  al30  be  seen  that  the 
oaloulations  describe  quantitatively  the  Turin  effect.  This  means  that  the 


method  can  simulate  polar  optical  scattering  as  well  as  i.ntarvalley 
scattering,  and  that  the  transition  from  polar  optical  scattering  (low 


snergy  region)  to  intervalley  scattering  (high  energy  region)  is 


accomplished  smoothly. 


Fig.  3-21  shows  the  calculated  electric  field  dependence  of  the  impact 
ionization  rate  in  GaAs  for  three  different  crystal  orientations.  We  have 
assumed  P  =  400.  The  shaded  region  indicates  the  range  covered  by  the 
experimental  data  (Section  3-2).  The  agreement  is  fair,  considering  the 
■uncertainty  in  the  scattering  rate  at  higher  energies.  The  inclusion  of 
upper  bands  i3  expected  to  increase  the  calculated  ionization  rats 
slightly,  and  therefore  to  improve  the  fit  to  the  experimental  data.  Mote, 
however,  that  the  calculation  shows  within  statistical  fluctuations  (~  20 
V)  no  orientation  dependece  for  the  ioinzation  rate.  This  contradicts  ths 
experimental  data  by  Pearsall  et  al .[39 ] ( ? ig ■  3-2).  Another  way  to 
calculate  the  orientation  dependence  is  by  rotating  the  electric  field 
direction  from  one  axi3  to  another.  The  result  is  shown  in  Fig.  3-22  to: 


an  electric  field  of  TOO  x’//cm  when  the 


f iel d  is  rotated 


<in>  to 


Electron  Ionization  Rate  (cm1) 


Fio-.  3.21.  Calculated  impact  ionization  rate  ot  an  electron  in  GaAs  with 
an  electric  field  in  three  crystallographic  directions  as  a 
function  of  a  reciprocal  field.  The  shaded  region  indicates 
the  range  of  available  experimental  data  (Fig.  3.1).  Within 
statistical  error,  we  do  not  see  any  orientation  dependence. 
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<1  1 1 >  direction.  Again  we  do  not  3ee  any  orientation  dependence  within  the 
statistical  error.  A  solution  including  this  rotation  pf  the  electric 
field  can  be  obtained  only  by  the  Monte  Carlo  method. 

The  effect  of  changing  the  value  of  P  has  been  al30  examined.  The 
result  is  shown  in  Pig.  3*23  for  the  case  of  ?  =  130,  400,  and  1 500  (the 
different  definition  from  Chwang[79]  has  been  used.  ?  =  100  in  our  case 
corresponds  to  P  =  50  in  their  paper) .  All  the  cases  give  practically  the 
same  impact  ionization  rate.  This  confirms  the  result  by  Baraff^s]  and 
Chwang  that  the  ionization  rate  is  insensitive  to  the  ionization 
probability  as  long  as  it  is  much  larger  than  the  probability  of  phonon 
scattering . 


Due  to  the  nature  of  the  Monte  Carlo  principle,  the  calculation  of  the 

effective  threshold  energyflOS]  is  trivial.  The  effective  threshold  energy 

is  defined  as  the  energy  at  which  Impact  ionization  actually  takes  place. 

It  is  usually  slightly  larger  than  the  threshold  energy  because  of  the 

finite  cross  section  for  pair  production.  The  result  of  t.nis  Monte  Carlo 

simulation  is  shown  in  Pig.  3-24.  The  effective  chreshold  energy  is  nearly 

independent  of  the  electric  field  for  large  values  of  P,  and  a  slightly 

increasing  function  of  the  field  for  smaller  values  of  ?  in  agreement  with 

others^??, 109] •  Thi3  also  agrees  with  the  result^  36 J  that  on  the  average 

the  imoact  ionization  event  takes  place  at  3.6  ±.3.2  of  aT\  above  the 

op 

threshold  energy  O  is  the  mean  free  oath  for  optical  oho  non  scattering-). 

op  ... 


Por  a  better  understanding  of  how  the  electrons  acquire  the  high 
energies,  and  how  impact  ionization  is  actually  accomplished,  we  show  in 
Pigs.  3.25  and  3*25  the  variation  of  electron  energy  after  each  scattering 
event  for  electric  fields  of  503  kV /cm  and  130  k 7 'em  respectively.  In  the 
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Fig.  3.25.  Variation  of  electron  energy  after  each  scattering  event  fo 
an  electric  field  of  500  kV/ctn  obtained  with  a  Monte  Carlo 
simulation. 
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Variation  of  electron  energy  after  each  scatterin 
event  for  an  electric  field  of  100  kV/cm. 
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case  of  500  kV./cm,  the  electron  energy  stays  around  0.7  ~  1.2  eV  most  of 
the  time,  'out  the  electron  occasionally  escapes  phonon  scattering  and  moves 
up  to  higher  energies.  In  Fig.  3*25  we  can  see  4  ~  5  spikes  which  reach  to 
~  1.8  eV.  When  an  electron  reaches  2.0  eV,  it  causes  impact  ionization. 
We  can  think  of  these  electrons  as  the  "lucky  electrons"  in  Shockley's 
theory,  and  those  electrons  around  the  average  energy  as  the  diffusing  (in 
energy)  electrons  in  Wolff's  therory.  However,  as  seen  in  the  figure  this 
classification  is  not  very  distinct.  Fven  those  electrons  in  the  spikes 
suffer  several  scatterings  before  they  reach  the  peak  energies.  Our 
results,  therefore,  contain  Shockley's  and  Wolff's  notions  of  ionizing 
electrons  as  does  Baraff's  theory,  but  under  much  more  general  conditions. 
Using  Baraff's  word!_7o],  the  notion  of  ballistic  electrons  by  Shockley  and 
diffusing  electrons  of  Wolff  are  "complementary"  in  determining  the  impact 
ionization  rate.  This  is  due  to  the  fact  that  the  height  of  spikes  (in 
Fig.  3.25)  depends  on  the  average  energy  of  electrons.  It  is  also 
important  to  note  the  difference  between  our  ionizing  electrons  and 
Shockley's  "lucky"'  electrons.  Shockley's  "lucky"  electrons  start  from  zero 
energy,  escape  the  phonon  scattering  completely,  and  impact  ionize.  The 
ionizing  electrons  of  our  result  3 tart  at  the  average  energy  and  reach 
ionization  threshold  after  a  few  scattering  eve>nt3.  This  explains  why 
Shockley's  theory  badly  underestimates  the  ionization  rate.  The  actual 
cause  of  impact  ionization  is  exactly  what  Shockley  neglected,  i.e.,  those 
electrons  which  suffer  phonon  scattering  at  intermediate  energies  and 
continue  upwards  in  energy. 
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3  -3  Summary 

A  Monte  Carlo  simulation  of  high  field  transport  in  GaAs  employing  a 
realistic  band  structure  has  been  described.  The  method  has  been  used  to 
study  the  impact  ionization  mechanism  in  Ga A3 .  The  band  structure  of  GaAs 
has  been  calculated  using  the  empirical  pseudopotential  method.  Partly  due 
to  the  lack  of  information  and  partly  for  simplicity,  we  have  made 
simplifying  assumptions  on  the  phonon  scattering  rates,  the  ionization 
threshold  energy,  and  the  ionization  probability.  This,  however,  is  not  an 
inherent  limitation  of  the  method.  Unlike  previous  theories  of  impact 
ionization,  the  method  requires,  in  principle,  no  adjustable  parameters  as 
long  as  the  band  structure  and  the  scattering  mechanism  are  known.  The 
method  has  provided  new  results  and  increased  the  'understanding  of  high 
field  transport  and  impact  ionization  in  GaA3.  The  calculated  drift 
velocity,  the  mean  free  path,  and  the  impact  ionization  rate  are  in  fair 
agreement  'with  the  experimental  data.  The  inclusion  of  the  higher 
conduction  bands  is  expected  to  further  improve  the  fit.  We  do  not  expect, 
however,  to  obtain  the  anisotropy  measured  by  Pearsall  st  al.[S9j.  In  our 
opinion  this  anisotropy  i3  not  a  consequence  of  the  band  structure,  but  is 
rather  caused  by  crystal  defects  or  other  effects  not  well  -understood.  It 
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phonon  scattering  rate.  The  effective  thereshold  energy  has  also  been 
calculated  and  discussed  briefly. 


Based  on  the  results  of  the  simulation,  a  general  discussion  of  impact 
ionzation  ha3  been  given.  We  find  that  typically  electrons  stay  around  an 
average  energy  and  experience  a  large  number  of  phonon  scatterings. 
Occasionally  electrons  escape  phonon  scattering  and  move  up  to  higher 
energy.  Some  reach  ionization  threshold  after  a  few  scattering  events. 
This  feature  is  seen  in  Fig.  3-2?.  It  can  be  considered  as  a  combination 
of  Wolff's  and  Shockley's  notion  of  ionizing  electrons,  but  the  distinction 
i3  rather  vague.  The  reason  for  the  success  of  3araff's  theory  is  that  his 
theory  also  contains  this  feature.  However,  because  of  his  formulation 
using  distribution  functions,  the  physical  picture  is  not  as  clear  as  in 
our  results.  Moreover,  our  method  includes  realistic  scattering  mechanisms 
and  band  structure  for  the  material  considered. 

Unlike  previous  theories  of  impact  ionization,  the  present  method  can 
in  principle  be  applied  to  any  semiconductor.  The  method  can  be  used  for 
both  polar  and  nonpolar  materials.  This  is  obvious  from  the  successful 
simulation  of  the  Ounn  effect,  which  contains  the  transition  from  polar 
optical  scattering  to  intervalley  scattering.  The  calculation  of  hole 
initiated  ionization  rates  should  also  be  possible,  although  presently  our 
understanding  for  hole  transport  is  not  as  deep  as  for  electron  transport. 
This  implies  a  tremendous  importance  for  understanding  the  operation  of 
photodetectors[3] ,  since  their  performance  depends  on  the  ratio  of  electron 
and  hole  initiated  ionization  rates.  It  should  be  understood  that  the 
method  is  quite  versatile  in  its  application.  A  transient  Monte  Carlo 
method  including  the  band  structure  may  be  used  to  investigate  the 
orientation  dependence  of  the  avalanche  response  time1"!  10  !.  The  inclusion 
of  the  position  dependence  should  enable  us  to  study  the  effect  of  the 
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There  are  some  improvements  which  should  be  investigated  in  the 
future.  First  of  all,  the  upper  conduction  bands,  at  least  the  second 
conduction  band,  should  be  included  in  the  simulation-  This  can  be  done 
rather  easily  if  we  know  (or  calculate)  the  scattering  rate  between 
different  bands.  The  inclusion  of  tunneling  is  also  possible.  Next,  finer 
me3h  points  should  be  used  in  the  calculation  of  the  band  structure  to 
increase  the  numerical  accuracy.  Ideally  the  maximum  energy  separation 
between  any  two  points  should  be  much  smaller  than  the  phonon  energy. 
Improvement  on  the  scattering  rates  also  seems  important.  The  scattering 
rates  should  be  assigned  at  each  k  point  including  the  overlap  integrals 
which  are  available  from  the  band  structure  calculations. 


124 


4.  CONCLUSION'S 

Two  theoretical  aspects  of  high  field  transport  in  III-V 
semiconductors  have  been  studied.  One  occurs  in  real  3pace  (real- space 
electron  transfer),  while  the  other  is  related  to  momentum  space  (band 
structure).  To  understand  and  exploit  high  field  transport,  we  have  to 
examine  it  in  both  real  and  momentum  spaces.  Just  as  the  band  structure 
(inhomogeneity  in  momentum  space)  makes  possible  the  transfer  of  electrons, 
inhomogeneities  in  real  space  (heterostructures)  lead  to  interesting 
transport  phenomena.  As  for  our  simulation  of  high  field  transport 
utilizing  a  realistic  band  structure,  we  feel  this  study  is  still  in  a 
preliminary  stage.  Our  ultimate  goal  is  to  understand  what  an  electron 
•undergoes  at  each  point  in  momentum  space  at  high  energies,  including 
details  of  the  scattering  mechanisms. 

I  would  like  to  close  this  thesis  work  by  quoting  a  statement  by 
3.  M.  Conwell[26]  in  1367,  "Our  present  knowledge  of  high-field  transport 
in  particular  materials  is  limited  to  regions  close  to  the  band  edges.  In 
the  near  future  we  should  be  able  to  penetrate  further  into  the  bands, 
proceeding  hand  in  hand  with  the  detailed  studies  of  band  structure  now 
becoming  available,  and  perhaps  even  contributing  to  these  studies.  Such 
advances  should  make  possible  a  more  detailed  unders tending  of  avalanche 
and  perhaps  ultimately  of  the  problem  that  first  stimulated  interet  'n 
these  studies,  dielectric  breakdown."  A  step  toward  these  goals  has  been 


made  in  thi3  thesis  work. 
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APPENDIX  1 

Total  Scattering  Rates  and  Angular  Dependences 
of  Various  Scattering  Mechanisms 

In  this  appendix  formulas  for  the  total  scattering  rates  and  the 
angular  dependences  of  various  scattering  mechanisms  for  charge  carriers 

— V 

with  wave  vector  k  are  summarized.  These  equations  are  needed  for  the 
Monte  Carlo  calculations  described  in  the  text. 

In  general  the  transition  rate  from  wave  vector  k  to  k'  is  given 
by  the  golden  rule: 

2 

S(k,k')  =  f1  |Hfi|  o(Ef-E.)  (Al.l) 


where  is  the  matrix  element  of  the  perturbing  potential  H  between  the 

initial  and  final  states.  The  angular  dependence  of  S(k,k')  gives  direct¬ 
ly  the  angular  dependence  of  the  scattering.  The  total  scattering  rate 
out  of  the  state  k  is  obtained  by  summing  over  all  the  final  state  k' . 
Therefore , 

00  TT  2tt 

-  =  — “  j  dk’  [  dp  da  S(k,k')  k'2  sin3  (A1.2) 

(2tt)  j  o  o  0 


where  V  is  the  volume  of  the  crystal,  and  3  is  the  angle  between  k  and  k' . 

RTierever  possible,  the  effect  of  non-parabolic  bands  and  the 
admixture  of  p-type  valence  band  wave  functions  in  the  T  valley  will  be 
included  in  the  formalism.  All  the  figures  include  non-parabolicitv ,  but 
not  the  admixture  of  p-type  wave  functions.  Non-parabolicity  is  defined 
by  the  relationship 


E(l+aE)  • 


(A1.3) 


Polar  Optical  Scattering  [33] 


The  transition  rate  for  polar  optical  scattering  is  given  by 


»  (  N  6  [E(k' )-E(k)-hw  ]  (absorption) 

S  (k,k’>  =  Q  B  •)  °  °  (A1.4) 

'  ’  '  ( (N  +1)6 [E (k ’ )-E(k)+hw  ]  (emission) 

o  o 


where 


2ire^h(a  f  ^ 

B  (k,k-)  =  - ++  2  7~  -  T~  G^,k’) 

4t re  v|  k-k’  |  ^  eo  £oo  • 

o 


(A1.5) 


In  these  expressions  is  the  optical  phonon  frequency  and  is  the 
thermodynamic  average  value  of  the  optical  phonon  number,  which  is  in 
equilibrium: 


o  exp(fu»  /k  T)-l 

O  D 


G(k,k’) 


1  f  ->■  *  ,-v 

-z  Z  |  dr  u  ->■  (r)u  ->.( 

y  5  p  1  J  y 1 k 1  yk 


(r)u  -j-(r) 


(A1.6) 


(A1.7) 


is  the  overlap  integral  between  the  periodic  parts  of  the  Bloch  functions 
at  k  and  k’  summed  over  the  doubly  degenerate  final  spin  states  and  aver¬ 
aged  over  the  initial  spin  states.  It  is  given  by 


where 


G (k, k’ )  =  (akak,  +  cRck,  cos  B)‘ 


l+aE(k)  j 
l+2aE(k)  j 


(A1.8) 


aE(k) 
l+2aE (k)  ! 


(A1.9) 


(A1.10) 


For  s-type  wave  functions  we  let  a  -*■  0  and  G(k,k')  =  1. 
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The  total  scattering  rate  i=  calculated 


where 


and 


em  "'ll) 


oil  1  j  (l+2ctE 1 ) 


4ms  /2  h 
o 


y"(e) 


F  (E,E') X 


\; 

\  *  o 


(Al.ll) 


E'  = 


F  (E,E* ) 
o 


E  +  ha>o  (absorption) 

j  E  -  hw  (emission) 

1  o 


C"1  [A&n  [Y,  (E)  1  Yr-E’ -[  +  B 
'Y^(E)  -  Y  ^ ( E ' ) ‘ 


(No+1) 


(A1.12) 


(Al .13) 


A  =  l  2(l+aE)(l+aE') +a(Y(E) +y(E,)>]' 


(Al . 14) 


B  »  -2aY^(E) y^(E ' )  [  4(l+aE) (1+aE*) 


+  ci{  y  (E)  +  Y  (E 1 )  }] 


C  =  4  (1+aE)  (l+aE ' )  (I+2aE)  (l+2ctE ' )  . 


(Al .15) 
(Al .16) 


To  recover  the  result  for  non-parabolic  bands  without  taking  account  of 
the  admixture  of  p-functions  we  let  a  -*■  0  in  A,  B  and  C,  then  A  =  C  =  4 
and  B  =  0 .  Quite  often  an  "effective  field  strength"  defined  as 


7% 

m  ehto 


Ll  _  a 


O  .  .2  Is  £ 

4me  n  >■  »  o 
o 


(Al . 17) 


is  used  in  the  total  scattering  rate.  It  has  the  value  of  Eq  =  5.95  kV/cm 
for  the  F  valley  of  GaAs . 

From  Eqs .  (A1.5)  and  (A1.8)  we  obtain  the  angular  dependence  of 


the  scattering 
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P(3)d6 


A  /  1+aE  A/l+aE' 

A  /  aE 

it  “E'  „ 

2 

_V  l+2ctE  V  l+2aE ' 

+  V l+2aE 

V l+2aE '  c  8_ 

sin  3  dg.  (A1.18) 


[Y(E)+y(E')  -  2Y"(E)y2(E')  cosS] 


To  realize  this  probability  distribution  using  uniformly  distributed  random 
numbers,  von  Neumann's  rejection  method  can  be  used  [33,38  ]  .  For 
parabolic  bands  with  a  =0,  we  have 


P(B)dS  ~ 


sin  3  d{3 


[  E+E'  -  2/EE1"  cos  3  ] 


(A1.19) 


In  this  case  using  random  numbers,  r,  equally  distributed  between  0  and  1 
gives 


S 


r  = 


P(3)dB 


(A1.20) 


or 


where 


cos  B  =  [(1+f)  -  (l+2f)  ]  /  f 


f  =  2/EEr  /  (/E  -  /!/) 2  . 


(A1.21) 


(A1.22) 


When  B  has  been  determined,  a  further  random  number  determines  the  azi¬ 
muthal  angle  between  0  and  2tt,  and  then  k'is  completely  specified. 

In  actual  Monte  Carlo  calculations,  however,  it  is  sometimes 
more  convenient  to  specify  the  phonon  vector  q  =  (for  the  case  of 

phonon  absorption)  by  the  vector  length  q  =  jq|  and  the  angle  9  between 
k  and  q  (Fig.  A.l).  Making  use  of 


q2  =  k2  +  k'2  -  2kk'  cos  3 


(A1.23) 


and  Eq.  (A1.21)  with 
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2/r(E)r(E') 

[/7(e7  -  /H^)]2 


(A1.24) 


appropriate  for  s-type  wave  functions,  it  is  possible  to  show  that 


q 


/y(E')  -  /y(E) 


y^r  (E 1 )  +  /r(E) 
/r(E')  -  /r(E) 


Once  q  has  been  determined,  cos  8  is  calculated  from 


k 


,2 


2 

q 


+  2qk  cos  0 


(A1.25) 


(A1.26) 


with  the  result 


~~  *  hm  [  1  +  a(2E  +h(i)  )]  +  q2 

ft2  o  —  o  M 

cos  9  =  - - — - - -  (Al.27) 

2q  •  ^  2m*y(E) 

where  upper  and  lower  signs  correspond  to  absorption  and  emission  of  a 
phonon  respectively.  Using  random  acimuthal  angles  relative  to  k  as 
before,  k'  again  can  be  specified. 

In  Fig.  A. 2  the  energy  dependence  of  the  polar  optical  scat¬ 
tering  rate  in  the  T  valley  of  GaAs  at  300°K  is  shown.  Fig.  A. 3  shows  the 
angular  dependence  of  the  scattering  probability  for  electron  energies  of 


0.1  eV  and  0.4  eV . 


Total 


Energy  dependence  of  the  polar  optical  scattering  rate  for  the 
valley  of  GaAs  (300  K) . 


Scattering  Probability  Density  (arbitary  units) 


Scattering  Angle  (radian) 

0  0.17T  0.27 t  0.37 r 


Fig.  A. 3.  Angular  distribution  of  the  polar  optical  scattering  probability 
for  electrons  with  energies  of  0.1  or  0.4  eV. 
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2.  Non-equivalent  Intervalley  Scattering  [33] 

We  distinguish  two  types  of  valleys  by  suffices  i  and  j.  The 
transition  rate  from  in  valley  i  to  "S'  in  valley  j  is 


S..(k,k') 


it 


N.  .6  [E.  (k')  -  E.  (k)  +  A  .  -A.  -flu..] 
ij  J  i  J  i 

(N.  .+1)5  [E.  (kf )  -  E. (k) + A .  -  A  .  +h  to  . 
ij  J  i  J  i  1J 


(A1 . 28  ) 


where 


B .  .  ( k ,  k ' )  -  Z.  -  h  -  D2.  G.  .  (k,k') 
i J  J  2pu  V  ij  ij 


(A1.29) 


and  a  is  the  intervalley  phonon  frequency,  is  the  number  of  valleys 
of  type  j,  D  is  the  intervalley  deformation  potential,  p  is  the  density 
and  j  is  the  intervalley  phonon  occupation  number  given  by  the  Bose 
distribution : 


ij  exp(fuo  .  . /k  T)-l 

ij  B 


(A1.30) 


The  energy  is  measured  from  the  minimum  of  the  valley  which  is  denoted  by 
A^  and  A^  for  valleys  i  and  j  respectively.  The  function  G„(k,k')  is 
given  by  an  equation  analogous  to  Eq.  (A1.8)  but  here  we  ignore  the  ad- 

-V  -> 

mixture  of  p-functions  and  let  G(k,k  )  =1.  The  total  scattering  rate  is 


,  v  3-, 

.  Z.DT.im.)  " 

1  =  .1  U  .1 

*  i  j  /2  to  h  ^  (  h'j  .  . ) 

ij 


(E')(l+2a.E') X 


(  (absorption) 

'  (N. .+1)  (emission) 
s  ij 


(A1.31) 
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where 


E.-A.+A.  +  hu>..  (absorption) 
i  3  i  ij 

E.-A.+A.-  fuo..  (emission) 
i  3  i  13 


(A1.32) 


All  the  intervalley  scattering  processes  including  equivalent 
intervalley  scattering  (explained  in  the  next  section)  are  completely 
randomizing,  i.e.  there  is  no  dependence  on  the  scattering  direction. 
Therefore  using  random  numbers,  r,  we  have 


p 

1 

=  —  sm  Sc 


(A1.33) 


cos  3  =  1  -  2r  . 


(A1.34) 


Fig.  A. 4  shows  the  energy  dependence  of  the  intervalley  scat¬ 
tering  rates  from  the  r  valley  to  the  L  and  X  valleys  of  GaAs  at  300°K. 


3.  Equivalent  Intervalley  Scattering  [33] 
The  transition  rate  is  given  by 


Se(k,k') 


-  f- ■><«-> 


N  6  [E(k')  -  E(k)  -hu  ] 
e  e 

(N  +1)6 [E(k')  -  E(k)  +  fh»>  ] 
e  e 


(A1.35) 


where 


Be(k,k’) 


(Z  -  1)  7  D  2 

e  2pw  V  e 
e 


(A1.36) 


and  w  is  the  intervalley  phonon  frequency,  D  is  the  deformation  potential, 
e  e 


Ga As  300  °K 
Intervailey  Scattering 


-  Total 

- F— X  Emission 

. T+L  Emission 

- f-^x  Absorption 

- f-^L  Absorption 


0  0.2  0.4  0.6  0.8  1.0  1.2 

Energy  (eV)  ,a_. 

Fig.  A. 4.  Energy  dependence  o£  the  intervailey  scattering  rate  in  the 
valley  of  GaAs  (300  K) .  The  L  and  X  valleys  are  0.33  and  0 
eV  awav  from  the  f  band  edge  respectively. 


and  Zg  is  the  number  of  equivalent  valleys.  Ng  is  the  thermal  equilibrium 
number  of  phonons  per  mode  given  by 


Ne  exp(hu  /k  T)-l 

e  d 


(A1.37) 


The  total  scattering  rate  is  calculated  as 


*  -72  2 

(m  >  De  V2  .  e 

(Z_-l)  — - o - —  Y  (E’)(l+2aE') 


N 


/Jirph  (haig) 


(A1.38) 


(N  +1) 


V.  '  e 


where 


E(k)  +  bcu  (absorption) 

.  i  _  )  fi 


(A1.39) 


1 E (k)  -  hrn  (emission) 
e 


In  Fig.  A. 5  the  energy  dependence  of  the  equivalent  and  non¬ 
equivalent  (L  valley  to  T  valley)  scattering  rates  in  the  L  valley  of  GaAs 
at  300°K  is  shown. 


4.  Acoustic  Phonon  Scattering  [33] 

We  assume  that  the  energy  of  the  acoustic  phonon  is  negligible 
compared  with  k^T,  and  neglect  the  change  in  the  electron  energy  during  the 
scattering.  Including  both  absorption  and  emission,  we  have 


S(£,k’)  =  ^  B  a  (£,£’)  N  6[E(k’)-E(k)] 

a  n  a  a 


(A1.40) 


where 


he 


2 

a 


!  k-k  ’  I 


G(k.k’) 


B  (k,£’) 


2 


T  «  —  t  r 


(A1.41) 


and 
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(■I 


hs  k-k' 


ft s  k-k' 


(equipartition) 


(A1.42) 


G(k,k!)  is  given  by  Eq.  (A1.8).  In  these  expressions  p  is  the  density, 
s  is  the  velocity  of  sound,  and  is  the  deformation  potential. 

The  total  scattering  rate  is  given  by 


•k  -fyo  2 

(2m  )  '  LT  S  2  , 

- -  y^(E)  (l+2aE)  F  (E) 

2^os2h4 


(A1.43) 


where 


F  (E)  = 


(1+aE) 2  +  ~(aE)2 
(l+2aE) 2 


(A1.44) 


To  recover  the  result  for  a  non-parabolic  band  alone  we  set  a  =  0  in  F  (E) . 

a 

The  angular  dependence  is  obtained  from  Eqs.  (A1.18),  (A1.38), 
(A1.39)  and  (A1.40)  with  the  result 


P(6)d8  ~  [ (1+aE)  +  aE  cos  3]  sin  8  dg 


With  random  numbers,  r,  the  angle  8  can  be  determined  by 


cos6  =  [ (l+2aE) 3 (l-r)+r]  /3  -  (14ctE) 

aE 


(A1.45) 


(A1.46) 


When  a  =  0,  the  scattering  is  completely  randomizing. 

In  Figs.  A. 6  and  A. 7  the  energy  dependence  of  the  scattering  rate 
and  the  angular  dependence  are  shown  for  the  T  valley  of  GaAs  at  300°K.  As 
seen  from  Fig.  A. 6  the  acoustic  phonon  scattering  rate  is  negligible  compared 
with  polar  optical  scattering  rate  and  intervalley  scattering  rate  in  polar 


Scattering  Rate  (sec-1) 


GaAs  300 °K 

Acoustic  Phonon  Scattering 


Fig.  A.  7.  Angular  distribution  of  scattering  by  acoustic  phonons  for 
electrons  with  energies  of  0.1  and  0.4  eV. 
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materials . 

5.  Intravalley  Non-polar  Optical  Scattering  [33] 

In  non-polar  materials  the  hamiltonian  matrix  element  of  non-polar 
optical  scattering  vanishes  for  s-like  band  minima.  This  is  also  true  for 
the  III-V  semiconductors  with  their  s-like  conduction  bands.  Consequently, 
this  process  needs  only  be  considered  at  high  energies.  At  these  energies 
it  is  probable  that  intervalley  scattering  dominates  in  most  of  the  III-V 
semiconductors.  Therefore,  non-polar  optical  scattering  is  commonly  ne¬ 
glected  . 

The  formula  for  non-polar  optical  scattering  is  very  similar  to 
the  one  for  intervalley  scattering.  The  total  scattering  rate  is  given  by 


_1_ 

t 

o 


where 


and 


2  *  V? 

D n  On  )  “  l/9 

y  (E ' )  (l+2ctE  ' )  X 


/Tuft  p  (hrn  ) 
o 


(v1} 


!E+h<u  (absorption) 

E  -  hm  (emission) 

^  o 


N 

o 


_ 1 _ 

exp (hrn  /k  T)  -  1 

O  D 


(A1.47) 


(A1.48) 


(A1.49) 


Dq  is  the  optical  deformation  potential. 

6.  Piezoelectric  Scattering  [27] 

In  cubic  crystals  lacking  inversion  symmetry,  the  piezoelectric 
stress  tensor  is  nonvanishing  and  electrons  are  scattered  by  the  piezo¬ 
electric  fields.  However,  piezoelectric  scattering  is  important  only  at 


low  temperatures  in  high-purity  materials.  Taking  into  account  screening, 
the  transition  rate  can  be  written  as 


2 it  fePz]  kBT 


S(k.k')  =  1 r  rrr1 


G(k,k' ) 


h  l££oi  2VC*  !k-k'|  +V2 


6  (E-  E’) 


(A1.50) 


where  P z  is  the  piezoelectric  constant,  =  ps  is  the  longitudinal  elas¬ 
tic  constant,  and  is  the  Debye  length  given  by 


ee  k_T 
o  B 

2 

e  n 


(A1.51) 


To  calculate  the  total  scattering  rate,  we  let  G(k,k')  =1  and  integrate 
over  the  final  state  k' .  Including  both  phonon  absorption  and  emission, 
we  get 


2  *  2 
1  a*  fePzl  kBX  1  +  2aE  .  f.  ,  S"  h>  ‘ 

pz  ^  o>  l  /2m\(E)  L  h 


(A1.52) 


From  Eq.  (A1.50)  the  angular  dependence  is  given  by 


P(B)d8  ~ 


sing  d  S 
h2 

[  y(E)  +  — *—5 y  (E)cosg  ] 

4m  Ld2 


(A1.53) 


Using  random  numbers,  r,  we  can  realize  this  distribution  as 


cosB  =  1  + 


4m  Ld2  y(E) 


1  -  1  + 


\  2  -1  r 

-Y~  y(e) 

1  _ 


(A1.54) 


If  one  wishes  to  determine  the  phonon  wave  vector  q  =  k'-k  first, 
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it  is  given  bv 


l  U 


L  v  u 


8m  Lq  y  (E) 


(A1.55) 


8  f 


and  cos#  is  calculated  as  in  Eq.  (Al.27) 


7 .  Ionized  Impurity  Scattering  [  35  ] 


The  transition  rate  is  given  by 


S(k.k')  = 


3  2  2 

32ir  (8e  )  N 


G(k,k’)_ 


0  ?  ? 

(4ir  ee  )  hV  (  jk-k  'I  “  +  \) 

o  1 


6(E- E'  ) 


(A1.56) 


where  is  the  impurity  concentration  and  \-  1/L^  is  the  inverse  screening 
distance.  G(k,k’)  is  given  by  Eq.  (A1.8).  The  total  scattering  rate  is 
calculated  as 


2n (Ze2) 2M, 


lim  (4 irzt  )  m 

0 


_ 1  +  2  a  E 

E3/2  (1  +  q  E ) 3 7 2 


:  [  1  +  2c  2 (X/ 2k)2]  2  o  2  ^  2 

x  - f - 7  +  ckMl  +  2ckZ%k)2]ln  t 

I  4(k/2k)-[l+(7./2k)-] 


O /2k) ~ 

1  +  ('‘./2k)  2 


where  c^  is  given  by  Eq.  (A1.10).  In  most  cases  of  interest 

2  2 

(2k/X)“  »  c^  .  Therefore,  we  get 


:II  (Ze2\2 


'  im  32  tfV2m 


(A1.57) 


_ 1 _ 1  +  2  qE 

(X/2k) 2  [1  +  (X / 2k)  2]  E3/2(l  +  qE)  3/2 


(A1.5S) 


For  the  angular  dependence  of  scattering  we  let  G(k,k)  =  1, 


fhen,  from  Eq.  (A1.56),  we  have 


By  using  random  numbers,  r,  we  get 


cos|3 


2(l-r) 

1  +  r  (4k2/X2) 


In  terms  of 


this  can  be  rewritten  as 


2  k  V'  1  -  r 
q  =  — — . — 

\Jl  +  r  (4k2/X2) 


and  cosd  is  given  by 


(A1.60) 


(A1.61) 


cosG  =  -  q/2k.  (A1.62) 

Fig.  A.-8  shows  the  energy  dependence  of  the  ionized  impurity 

17  -3 

scattering  rate  for  uncompensated  GaAs  doped  with  NQ  =  1  x  10  cm  and 

18  -3  18  —3 

for  compensated  GaAs  doped  with  =  1.1  x  10  cm  and  =  1  x  10  cm 

17  -3 

which  gives  an  electron  density  of  1  x  10  cm  .  In  uncompensated  material 

the  scattering  rate  depends  very  little  on  the  impurity  doping  density  due 

to  screening.  Fig.  A. 9  shows  the  angular  dependence  of  ionized  impurity 

scattering  for  electron  energies  of  0.1  eV  and  0 . 3  eV  and  two  different 

16  ~3  17  -3 

electron  densities,  1  x  10  cm  and  1  x  10  cm  As  can  be  seen  from 
this  figure  ionized  impurity  scattering  gives  small  scattering  angles. 

8.  Alloy  Scattering  [73] 

Here  we  adopt  the  model  developed  by  Harrison  and  Hauser.  They 
assume  the  crystal  to  have  a  uniform  background  potential  with  square  wel] < 
of  depth  AE  and  extent  r  <  r^  at  random  sites  associated  with  one  of  the 
alloying  consistuents .  They  arbitarily  take  rQ  to  be  the  nearest-neighbor 


Scattering  Rate  (sec  L) 


Ga  As  300°K 

Ionized  Impurity  Scattering 


Compensated 

N0  =l.lxl018crrr3 
NA-lxl018cm'3 
n  =lxl017cm'3 


Uncompensated 
N0=lxl017cm"3 
n  =lxl017cm'3 


0  0.2  0.4  0.6 

Eneray  (eV) 

Fig.  A. 8.  Energy  dependence  or  ionized  impurity  scattering  in  uncompensated 
and  compensated  GaAs  (300  K) . 


Scattering  Probability  Density  (arbitrary  units) 


Fig.  A. 9.  Angular  distribution  of  scattering  probability  by  an  ionized 
impurity  for  electrons  with  energies  of  0.1  and  0.4  eV  and 

16  -.3 

electron  densities  of  1  x  10  cm  and  1  x  10x  cm"  . 
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separation,  or  -~V  3  a.  AE  is  taken  to  be  the  electron  affinity  difference 
between  two  end  binaries.  Assuming  the  perfectly  random  case  gives,  for 
the  transition  rate, 

2 


S  (k,k')  =  ^  .  -y^Tv  •  X(1-X)  (AE)2  5(E_E'  } 

L  a 


(A1.63) 


where  N  is  the  density  of  alloying  sites,  and  x  is  the  mole  fraction  of 

3. 

the  binary  AC  in  the  ternary  alloy  A^B^  ^C. 

Integrating  over  the  final  states  gives  the  total  scattering 


rate 


-v /  ’  2  *  3/2  1/2 

1  =  3  A/  2 if  .  x  (1  -  x)  (AE)“  (m  )  /_  (1  +  2sE)t  '  (E) 

Tal  8h4^ 

a  (A1.64) 

As  seen  from  Eq.  (A1.63)  alloy  scattering  is  completely  randomizing  in  the 
first  order.  Therefore,  cos3  is  determined  by  Eq.  (Al.34). 

Fig.  A. 10  shows  the  energy  dependence  of  alloy  scattering  rate 
for  Aln  0  Gan  Q  As  (x  =  0.2)  at  300°  K.  (m  =  0.0796  m  ,  AE  =  0.55  eV  and 

(J  •  Z  U .  o  O 

22  -3 

N  =  2.22  x  10  cm  have  been  used  for  Che  calculations.) 
a 


Energy  (eV) 


Energy  dependence  of  random  potential  alloy  scattering  rate  in 
the  valley  of  A1q  2Ga0  8As  K) •  The  model  is  by  Harrison 

and  Hauser  [73  ]  . 
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APPENDIX  2 

Material  Parameters  used  in  the  calculations  [34  ] 

1 .  GaAs 


Bulk  material  parameters 

Lattice-  constant 

5.642  X 

Density 

5.36  g/cm8 

Electron  affinity 

4.07  eV 

Piezoelectric  constant 

0.16  C/m2 

LO  phonon  energy 

0.03536  eV 

Longitudinal  sound  velocity 

5.24  x  105 

cm/ sec 

Optical  dielectric  constant 

10.92 

Static  dielectric  constant 

12.90 

Valley-dependent  material  parameters 

r(ooo) 

L (111) 

X(100) 

J. 

Effective  mass  (m  /m  ) 

o 

0.063 

0.222 

0.58 

Nonparabolicity  (eV 

0.610 

0.461 

0.204 

Energy  band  gap  relative 
to  valence  band  (eV) 

1.439 

(0) 

1.769 

(0.33) 

1.961 

(0.522) 

Acoustic  deformation 
potential  (eV) 

7.0 

'9.2 

9.7 

Optical  deformation 
potential  (eV/cm) 

0 

3  x  108 

0 

Optical  phonon  energy  (eV) 

— 

0.0343 

— 

Number  of  equivalent 
valleys 

1 

4 

3 

Intervalley  deformation 
potential  (eV/cm) 

r 

0 

9 

1  x  10s 

1  x  109 

L 

9 

1  x  10 

1  x  109 

5  x  108 

X 

9 

1  x  10 

5  x  108 

7  x  108 

Intervalley  phonon 
energy  (eV) 

r 

.0 

0.0278 

0.0299 

L 

0.0278 

0.0290 

0.0293 

X 

0.0299 

0.0293 

0.0299 

2.  A1  Ga..  As 
x  1-x 

Effective  mass  in  the  f  valley 

Energy  bandgap  (eV) 

Optical  dielectric  constant 

Static  dielectric  constant 

m  /m  *  0.063  +  0.083x 

0 

Er  =  1.439  +  1 . 247x 
g 

zm  =  10.92  -  2.42x 

e  =  12.9  -  2. Ox 
o 

LO  phonon  energy  (eV) 


Hu  =  0.03535  +  0 . 014 64x 
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APPENDIX  3 


Wave  Vector  Relations 


In  this  appendix  the  vector  algebra  necessary  to  determine  k' 

— ^ 

is  presented.  kf  can  be  determined  either  directly  by  specifying  the  polar 

— ►  — ►  — *  |  — ► 

coordinates  (k  ,  3,  a)  relative  to  k,  or  by  specifying  q  =  k  -  k.  We 
write  the  vector  k  as  (refer  to  Fig.  A. 11) 


/  \\  f  ksin#  1cos4>1  \ 
k  =  k  i  =  f  ksin#  sin$1  \ 

\  y  ■  \  1  1 

\kz  /  \  kcos#  x  / 


(A3 .1) 


1.  k1  specified  by  the  polar  coordinates  (k1,  3,  a)  relative  to  k 

The  length  of  k'  is  determined  by  energy  conservation.  We 

— V 

consider  an  operation  which  aligns  the  vector  k  with  the  z-axis.  This 
operation  can  be  written  as 


U  *  Y  Z 

“®1  “*1 


(A3. 2) 


where  Y  represents  the  rotation  around  the  y-axis  with  an  angle  8  ana 
6 

Z  represents  the  rotation  around  the  z-axis  with  an  angle  b.  in  matrix 

4> 

representations  these  rotations  can  be  written  as 


/cos 9  0  sin0  \ 

YQ  =  [  0  1  0 

\-sin0  0  cos 0  . 

I  cos$  -sin6  0  '• 

Z  j.  =  !  sinb  cos6  0 

$ 

0  0  1/ 


(A3. 3) 


(A3. 4) 


Now  suppose  the  angle  3  between  k  and  k’  and  the  azimuthal  angle,  a, 

— ^ 

are  given  as  in  the  Monte  Carlo  calculations.  Then  k*  is  given  by 


— *  _  1  -►  -+ 

k  ’  =  U  k  "  =  Z .  Y0  k  ” 

01  9l 


(A3. 5) 
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1 


«  r 


where 


'k'  sinScosa\ 
k'  '  =[  k’  sinSsina  j 
^k'  cos8  J 


li:  “ 


SB  r 


r 


Performing  the  matrix  multiplications,  we  obtain 

kx*  =  k'  (cost)  ^cos<J>^sinScosa  -  sin<{>^sin3sina  +  sin0^cos4>. 

k^'  =  k'  (cosfl  ^sin^singcosa  +  cos^singsina  +  sin#  sin^., 

k  '  =  k'  (-sin# , sinScosct  +  cos0,cosS). 
z  1  1 

->  — *  — ► 

2.  k  specified  by  q  =  k'  -  k 

— ► 

Given  the  length  of  q  and  the  polar  angles,  6  and  6, 

— V  — V 

to  k,  q  can  be  obtained  by  operating 


/  cosfl,cos4  -sm?,  sm0,cos6i  \ 
ii  l  l  i  \ 

r 

‘S,  =  '  cos$  l^in^q  cosb-^  sinfi^sinb-^ 


Z,  Y. 

1  li 


to 


\  -3in9l 


/  q  '  \ 

x  \ 


q’  - 


0  cos0 

qsin0cos£  \ 

\ 

qsin^sind  I 
qcos4> 


Therefore,  using  k'  =  k  +  q,  we  have 


k^’  =  q^'cos#  ^cos<j)^  -  q  'sin<j>^  +  (q^  ’+k) sin#  ^  cos<{> ^ 
k  '  =  q^ '  cosS  ^sin:f>  ^  +  q^'cosb^  +  (qz'+k)  sir.0^sini^ 
+  (qz'+k)cos31. 


(A3. 6) 

cosB) 

cosB) 

(A3. 7) 

relative 


(A3. 8) 


(A3. 9) 


(A3. 10) 
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